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■^j- ■ We propose a nonperturbative approach to nonabelian two-form gauge theory. We formulate the 

CN 1 theory on a lattice in terms of plaquette as fundamental dynamical variable, and assign U(N) Chan- 
O 

Paton colors at each boundary link. We show that, on hypercubic lattices, such colored plaquette 
variables constitute Yang-Baxter maps, where holonomy is characterized by certain dynamical 
deformation of quantum Yang-Baxter equations. Consistent dimensional reduction to Wilson's 
lattice gauge theory singles out unique compactness condition. We study a class of theories where 

c$ ' 

the compactness condition is solved by Lax pair ansatz. We find that, in naive classical contin- 
uum limit, these theories recover Lorentz invariance but have degrees of freedom that scales as 
A^ 2 at large iV. This implies that nontrivial quantum continuum limit must be sought for. We 
demonstrate that, after dimensional reduction, these theories are reduced to Wilson's lattice gauge 
theory. We also show that Wilson surfaces are well-defined physical observables without ordering 
ambiguity. Utilizing lattice strong coupling expansion, we compute partition function and correla- 
tion functions of the Wilson surfaces. We discover that, at large N limit, the character expansion 
coefficients exhibit large-order behavior growing faster than exponential, in striking contrast to 
Wilson's lattice gauge theory. This hints a hidden, weakly coupled theory dual to the proposed 
tensor gauge theory. We finally discuss relevance of our study to topological quantum order in 
strongly correlated systems. 



1 Pictures and Discussions 



"What is the use" thought Alice "of a book if it has 
no picture or conversation ?" — Louis Carroll 

An outstanding problem in theoretical physics is a constructive definition of p-form gauge theories, 
especially, nonabelian and self-interacting ones. Variants of p-form gauge theory arise in diverse 
contexts, ranging from string or M theories (Bill and higher-dimensional integrable systems 
to topological order and phases in strongly correlated systems [4] and to quantum error correction 
codes [[51 in quantum information sciences. Of particular interest is whether a nonabelian p-form 
gauge theory exists and, if so, what sort of self-interactions are allowed by the gauge invariance. 
To the problem posed, one's first guess is that the fundamental degrees of freedom are some sort 
of nonabelian extension of the abelian p-form gauge theory, but then the question is now more to 
"what types of nonabelian extensions can be endowed to abelian p-form gauge theory?" and to 
"what types of self-interaction are possible for a given nonabelian extension?". 

The Z2 Ising model provides the simplest situation of all. Consider the model defined on 
a d-dimensional hypercubic lattice. In dual formulation, the Ising model is mapped to a gauge 
theory, where the gauge potential is a p = (d — 2) -form and takes a value in G = Z2. As is firmly 
established, the Ising model does not exhibits any nontrivial renormalization group fixed point for 
d > 4. It implies that the dual p-form gauge theory in d > 4 ought to be free in the continuum 
limit, yielding no obvious self -interaction among the p-form gauge potentials!*]. 

In string and M theories, variety of the p-form gauge potentials is rich and complex. Foremost, 
all known string theories, whether supersymmetric or not, contain universally the Kalb-Ramond 
two-form potential B2 = jjB mn x m A x", and the fundamental string couples minimally to it [1 1 . Type 
I and II superstring theories contain in addition R-R(Ramond-Ramond) p-form potentials [2]. The 
D(p — l)-branes are the charged objects coupled minimally to these R-R potentials. Because of 
different chirality projection, type IIA superstring gives rise to p = 1,3, ■■■,9 only, while Type IIB 
superstring does so for p = 0, 2, • ■ ■ , 10 only. These NS-NS and R-R fields exhaust all possible p- 
form potentials permeating through the ten-dimensional bulk spacetime. To the extent understood 
so far, they are essentially abelian. 

There are, in addition, p-form potentials residing only on the worldvolume of string solitons. 
Consider first the D-branes, the objects minimally coupled to the R-R gauge potentials. At low- 
energy below the string scale, D-brane worldvolume dynamics is dominated by the lowest exci- 
tation of open strings whose both ends are attached to the D-brane. The excitation constitutes 
one-form (p = 1) potential A m (x)x of gauge group U(l) and free scalar fields. A novelty is that, 

*The Ising fixed point exhibits strong stability under renormalization group flow. For instance, even dense dilution 
over simplices of nonzero codimensions is unable to induce a flow away from the Ising fixed-point Q. This implies 
that the dual p-form gauge theory remains free even for randomly distributed coupling parameter. 
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when N parallel D-branes stack on top of one another, combinatorially, there are N x N possible 
open strings and the lowest excitation of them forms U(AT) matrix-valued 1-form potential |[3). 
Thus, for D-branes, N 2 variety of open strings constitutes the microscopic field and particle de- 
grees of freedom of the D-brane worldvolume. Consider next the NS5-brane, the magnetic dual to 
the fundamental string. Worldvolume dynamics of an NS5-brane in Type IIA string theory, equiv- 
alently, an M5-brane in M-theory is described at low-energy by the six-dimensional 9\£ = (2,0) 
superconformal theory lfT8l . and the theory is known to contain a selfdual two-form (p = 2) po- 
tential of gauge group U(l). Again, a novelty is that, when N parallel NS5-branes or M5-branes 
stack on top of one another, the microscopic degrees of freedom constitute tensionless strings that 
arise in M-theory from the variety of open M2-branes connecting all possible combinatoric pairs 
of the N M5-branes iflOl . Similar to the D-branes, one might anticipate that N 2 variety of open 
M2-branes connecting all possible pairs of M5-brane constitute the microscopic field and string 
degrees of freedom of the M5-brane worldvolume. However, in stark contrast, various consid- 
erations ranging thermodynamic free energy ifTTTl and gravitational anomaly cancellation lfT2l all 
indicate a peculiarity that the degrees of freedom are intrinsically quantum-mechanical and scales 
in N — > oo limit as N 3 , in stark contrast to N 2 behavior ITT31 observed for the D3-brane worldvolume 
dynamics described by the 2\£ = 4 super Yang-Mills theory. 

In this paper, we study a viable extension of nonabelian gauge invariance to tensor gauge field 
and put forward a nonperturbative approach for tensor gauge theory in d > 4 dimensions Q by 
putting the theory on a J-dimensional hypercubic lattice. The lattice formulation allows a tran- 
parent and concerete description for origin of nonabelian gauge symmetries and self-interactions 
among the p-form gauge fields By putting the theory on a lattice, we are sidestepping from 
other structures that can be endowed to the tensor gauge theory such as (extended) supersymmetry 
or (anti)self-duality. Though these structures are desirable for making contact with those arising 
in string theory, in this paper, we focus primarily on the issue of nonabelian gauge symmetry and 
self-interactions thereof. 

This paper is organized as follows. We begin in section 2 with etiology of Chan-Paton factors. 
We assign Chan-Paton factors to boundary links of an elementary plaquette so that it carries four 
'color' indices. We take these objects as fundamental dynamical variables and construct in section 
3 a nonabelian two-form tensor gauge theory defined on a d-dimensional hypercubic lattice. Ex- 
pressing the plaquette variable as (N 2 x N 2 ) matrices of U(A0 gauge group, we construct action 
for nonabelian tensor gauge theory. We then study possible compactness conditions. We show that 
consistent reduction to Wilson's lattice gauge theory lfT6ll by a dimensional reduction and unitarity 
or reflection-positivity singles out a unique choice of the condition. In section 4, we study ground 



^In this paper, we study exclusively two-form gauge theory — as will become evident in foregoing discussions, 
the construction is extendible to higher /5-form gauge theories straightforwardly. 

^There has been in the past occasional attempt for constructing nonabelian p-form gauge theories. See lfT4ll and 
also 031. 
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state of the lattice theory. We show that ground-state configurations are provided by vanishing 
holonomy of the nonabelian tensor fields. Remarkably, these configurations are specified by the 
solution of so-called dynamical quantum Yang-Baxter equations, pointing to an extremely rich 
structure of the ground-state wave function. Nontrivial holonomies are measured by a set of possi- 
ble deformations of the dynamical quantum Yang-Baxter equations and their solutions. In section 
6, we study Lax pair ansatz for solving the compactness condition by parametrizing the plaquette 
variables in terms of direct product of (N x N) matrices of U(A0 group. We then study continuum 
limit of the classical action, viz. naive continuum limit and point out that the gauge symmetry 
becomes abelian in this limit. In section 6, we show that the theory passes up an important consis- 
tency condition: upon lattice dimensional reduction, the theory is reduced to the Yang-Mills theory 
in the classical continuum limit. In section 7, utilizing the lattice strong-coupling expansion, we 
compute free energy and Wilson surface correlators. We find that in stark contrast to Polyakov- 
Wilson lattice gauge theory, the expansion series is not Borel summable even in the large N limit. 
We present an intuitive argument for the behavior and argue that the large-order behavior point to 
the existence of a weakly interacting, dual lattice theory. Section 8 is devoted to a summary of 
this paper and discussion on remaining issues. In Appendix A, we explain that a truncation of the 
internal degrees of freedom following the dimensional reduction leads the Wilson's lattice gauge 
theory at the lattice level the plaquette variables to the Wilson's link variables. Consistency with 
this procedure reduces the original four possibilities of the theory into two. In Appendix B, we 
analyze the rank-./V theory and see that it does not lead a physically meaningful continuum theory 
under the dimensional reduction. In Appendix C, we explain solutions of the compactness condi- 
tion of the rank-/V 2 (I) which are not expressed as a direct-product structure. Although the solutions 
have the (N 3 ) degrees of freedom, the gauge degrees of freedom reduce to (11(1))^, which does 
not lead to an interesting nonabelian tensor theory. Also, since the gauge degrees of freedom are 
not enough to eliminate modes of wrong-sign kinetic terms (ghosts), the solutions do not seem to 
lead physically meaningful continuum theory at least in the classical level. In Appendix D, we 
present some computation on a character expansion coefficient used in strong coupling expansion. 
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2 Etiology of Chan-Paton Factors 



We begin with etiology of the Chan-Paton factors [fTTl . Originally, the Chan-Paton factors were 
introduced as a prescription for introducing "colors" to open string amplitudes. In this section, 
we will adapt the notion on a lattice and interpret it as endowing a Chan-Paton bundle over a 
finite-dimensional vector space V*cp to parallel transport. 

2.1 elementary link variables 

Consider the Wegner-Wilson-Polyakov formulation of the lattice gauge theory lfT6ll on a J-dimensional 
hypercubic lattice. Taking gauge group G to be compact U{\), elementary degrees of freedom are 
associated with link variables, U(x,x + jj). This is a dynamical variable residing at each unit link 
(of length a) [x;fi] of direction ft, connecting the nearest neighbor sites x and x + jj. The link vari- 
able assigns a dynamical weight to parallel-transporting a charged particle (charge +1) long the 
unit link [x;fi\. The link variable U^(x) is parametrizable in terms of gauge potential A jU (x) that 
takes values on § : 

Ufj(x) = U(x,x+jj) = exp(i'oA (U (jc)), (2.1) 

By the exponential map defined so, the link variable satisfies the compactness condition: 

U- l (x,x + ft)=U(x + /j,x). (2.2) 

As the charged particle is parallel-transported, U^(x) creates the electric charge —1 at the site x 
(where a unit charge is depleted) and +1 at the site x + jj (where a unit charge is deposited). For 
an arbitrary path on the lattice, the corresponding parallel-transport is determined by multiplying 
link variables along the path. Therefore, the set of all possible link variables {U^(x)} constitutes 
microscopic degrees of freedom on the lattice. Of particular interest is the parallel-transport around 
a unit plaquette [x;//,vj: 

Up = Uix.x + ^Uix + fj.x + fi + ^Uix + fi + v.x + ^Uix + v.x). 

This is a gauge-invariant operator, and deviation of it from unity measures the holonomy around 
the unit plaquette v]. 

The construction is readily extendible by adjoining Chan-Paton factors. We promote the link 
variable to a matrix- valued one. It comes through a pair of the Chan-Paton factors G Vcp 
attached at the two ends of each unit link at x: 

U M (x)®\ij)= (fyto) (2-3) 

The ij= 1, • • • ,N indices label an orthogonal basis of the vector space V*cp, so equivalence rela- 
tions and hence the gauge group G would be associated with the rational map 

G : Vcp — > Vcp- 
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As such, the elementary link variables are matrix- valued in GL(Af, C). See FigJU 

The Chan-Paton etiology can be restated as follows. The sites and links are 0- and 1-simplices 
of the J-dimensional hypercubic lattice. The two sites at the end of a link variable are boundary 
0-simplices of a 1 -simplex. To define nonabelian lattice gauge theory, we are then assigning the 
Chan-Paton factors z',y, • • ■ G Vcp at each boundary 0-simplex of 1-simplices. The matrix-valued 
holonomy is measured around the boundary of each plaquette. For the unit plaquette located at x 
and oriented along [/iv] -directions, the holonomy is measured by 

(^W) y (Cv(,+A)). i (^+A+*)) K (f/-v(,+*)) <m 

It measures the gauge flux felt by a colored particle around the closed loop of 1-simplices, trans- 
muting the color through i — > j — > k — >■ £ — >■ m. In the context of open strings ending on D-branes, 
the colored particle is nothing but an endpoint of the open string. 









1 v J 


c 

i x 


3 > • A 

.V+U 







Figure 1 : The link corresponding to the variable (Ufj(x))jj. The Chan-Paton indices i, j are assigned 
as quantum numbers labelling the 0-simplices at the two ends of the 1-simplex defining the link 
variable. 

The elementary link variable (UJx))ij evolves dynamically and sweeps out a plaquette P (of 
area a 2 ) inside the J-dimensional hypercubic lattice. The Chan-Paton factor should be consistent 
with various closure relations. First, a product of two link variables should be isomorphic to 
another link variable: 




(2.4) 



Second, generalizing Eq.(??), unitarity of the link variable forces the compactness condition: 




(2.5) 
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This reduces the space of link variables from GL(N, C) to U(N). Third, the gauge invariance 
requires that 

( f /„ W ).^(v W ) a ( f /, W ) H (vt(, + ?)) ( . (2.6) 
is an equivalence relation. 

2.2 elementary ^-simplex variables 

The above consideration is readily generalizable to higher simplices. For link variables, crux of 
the idea to the Chan-Paton factor was that 'color' degrees of freedom is attached to the boundaries 
(O-simplices) of the link variables (1 -simplices). As the boundary consists of two elements, the 
resulting Chan-Paton gauge groups are simply matrix group GL(N). We shall now extend this 
notion to higher-dimensional simplices, viz. assign Chan-Paton indices to the (p — l)-simplicial 
boundaries of a p-simplex. 

Consider again rf-dimensional hypercubic lattice, and take p-simplices A p (2 < p < d — 1) 
and p-pairs of (p — 1) -simplices A p _i instead of 1 -simplices and one pair of O-simplices, respec- 
tively, and assign a dynamical variable U Pl ... Pp (x) on the /^-simplex located at x and oriented along 
\pl, ■ ■ ■ ,jjp\. The set of these p-simplex variables constitute the microscopic degrees of freedom 
of p-form lattice gauge theory. So, for the compact U(\) gauge theory, we may parametrize the 
p-simplex variable by a direct extension of Eq. (l2. II) as 

U n ... Pp (x) = exp (ia p A ltl ... Pp (x)) 

viz. an exponential map for the variable A pi ... Pp (x), and interpret it as the dynamical weight for 
parallel transporting a charged (p — 1) -dimensional brane by lassoing it around the p-simplex A p . 
Extending the interpretation, one can now assign the Chan-Paton factors along the boundary of the 

S -simplex variable, viz. the (p — l)-simplices forming the boundary of p-simplices: A p _i = dA p 
, Assigned this way, gauge transformation of the nonabelian p-simplex variable is defined as 

u lxv .. IXp (x)^a(x)-u n ... flp ( x )-a- l (x), 

where the multiplication is defined in terms of appropriate action on the p-tuple of the Chan-Paton 
vector space, Vcp <8> ■ ■ • <8> Vcp. The curvature of the p-form gauge potential is measured by the 
deviation of 

( n n ow*)) 1 

§More generally, one can assign (p — 1) -dimensional "boundary" field theory living on A p -\, whose excitations 
would carry a "continuous color" degrees of freedom in replacement of the finite-dimensional Chan-Paton indices. 
We thank E. Witten for suggesting this extension. 
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from unity, where again multiplication of variables is defined appropriately on the p-tubles of the 
Chan-Paton vector space, £g> p Vcp- Evidently, the holonomy and the curvature are valued in the 
(p + 1) -rubles of the Chan-Paton vector space, ® p+1 Vcp. 

2.3 Universality 

Having generalized as above, one immediately realizes that there is an enormous difference be- 
tween the link variables and the higher-dimensional extensions. In discretizing the space, we have 
tacitly chosen d-dimensional hypercubic lattice. We could have chosen an alternative lattice lead- 
ing to a different discretization. Then, among all p-simplex variables, the link variable is special 
since, irrespective of the lattice type, boundary of the link variable consists always of two points. 
This is essentially a statement of the 'universality' — detailed specification of the discretization 
should not matter for long-distance and continuum limit of the theory defined so. 

Thus, universally, a link variable with Chan-Paton indices would give rise to N 2 degrees of 
freedom. For higher-dimensional simplices, the number of boundary simplices depend explicitly 
on details of the latticization. For instance, consider two types of plaquettes: square or triangle 
ones. For each of them, we would be assigning four or three Chan-Paton indices. Then, in sharp 
contrast to link variables, a plaquette variable of 3 or 4 Chan-Paton indices may give to iV 3 or N 4 
degrees of freedom depending on the choice of latticization. Interestingly, the degerees of freedom 
is at least N 3 , as the minimal choice is lattice with triangular plaquettes. 

Remarkably, in the next section, we will find that compactness conditions, which are direct 
generalizations of the condition Eq. (l2.5l) . constrain it so severely that the plaquette variable actually 
carries o(N 2 ) degrees of freedom. 

3 Nonabelian Tensor Gauge Theory on a Lattice 

In this section, built upon the result of the previous section, we construct a nonabelian tensor gauge 
theory defined on a J-dimensional hypercubic lattice. 

3.1 lattice action 

We now utilize the observation we made concerning the Chan-Paton prescription to a plaquette 
variable U^ v (x), defined on an elementary plaquette located at site x along /j,v directions See 
Figl2l Heuristically, the 2-simplex is interpretable as the trajectory of a string on the plaquette 
[x 7 x + ju,x + /u + v 1 x + V,x], and the plaquette variable is the amplitude of parallel-transporting a 
"colored" string from two adjacent links along /},v directions to diagonally opposite links. There 
are two possible choices of the "move" but they are related by the discrete subgroup of 0(d) 
rotation group. The boundary of the 2-simplex consists of one-simplices: [x,x+fi] : [x + £i,x + fi + 
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v], [x + jj + V, x + v], [x + V,x]. We then attach Chan-Paton factors /, j,/c, I to each 1-simplex, taking 
values in the newly introduced 'color' space Vcoioj-. This 'color' dressed 2-simplices will be the 
basic building blocks of the nonabelian tensor gauge theory we are proposing. In the followings, 
we takeV co i or = U(AO. 

Consider J-dimensional hypercube lattice, whose sites are labelled as x = (x , ■ ■ - ,xr) and unit 
vectors are labelled as //, V,X, • • ■. The action we propose is a direct generalization of the Wegner- 
Wilson-Polyakov prescription: 



d 



Hensor = "P £ E ' Re [ U \- X '^ ~ N 
{x}/j,v,X=1 



(3.1) 



where P denotes inverse coupling parameter, and 11 [C^M] refers to the gauge-invariant action 
density for an elementary cube [jc;//vX], defined as 

N 

u[x-m] = £ (u^(x)j rki (u vX (x)) im (u^(x)^ 



(3.2) 



xfu, v (x + X)Y (u vX (x + ju)Y (U^(x + V)Y 

\ / rstu \ / jpsm \ / qtpk 

(3.3) 

= Tr(([/-C/.C/) comer (x)-([/-[/-C/) c t omer (x)). (3.4) 

The last expression is to emphasize the interpretation that the action density is a 'trace' over the 
string holonomy, viz. phases acquired through three adjacent plaquettes on a corner located at x and 
their hermitian conjugates. In Eq. (l3.1l) . the (') primed sum over fj, v, X runs over the J-dimensional 
lattice directions, where no pair of them are allowed to coincide. The variable {U^,{x))ijki is 
complex- valued, and lives on the plaquette encompassing the sites x, x +//, x +p + V, x+ V, x. The 
Chan-Paton indices z, j, /c, / are associated to the four links forming the boundary of the plaquette 
as depicted in Fig. [2l The complex conjugation V refers to simultaneous reversal of the plaquette 
orientation and the Chan-Paton colors at the boundaries. As such, we adopt the convention that the 
plaquette variables satisfy 



E^(x)J. = (U VfJ (x) (3.5) 

/ ; jkl \ / Ikji 

Notice that ordering of both the directional indices /j, v and the color indices i, j, /c, / are reversed 
by the complex conjugation. 

We also define a gauge transformation rule as 



-w)„ H ^ , E,( v 'W)»(«'^ + »)^( l ''-w), /w ( v ^ + ^ t )«K) t ) n 

= ((V^-VvlV-^-Vj)^)) , (3.6) 
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Figure 2: The plaquette corresponding to the variable ( U /M (x) ) for n ^ v. The Chan-Paton 
indices are assigned on each link (1 -simplex ) as i, j, k, I = 1, — ,iV. 

where V^(x) are gauge transformation functions valued in U(N) group, obeying V/j(x)V^(x) = In- 
It is easy to see that 11 [C^-^x)] in the action S is associated to a 3-simplex (cube) as in Fig. [3l and 
invariant under the gauge transformation Eq. (l3.6l) . 



3.2 compactness conditions 

In Wegner-Wilson-Polyakov lattice gauge theory, the link variables (U^(x))ij are assumed com- 
pact, as defined via the conditions 




(3.7) 



It means that (U^{x))ij are represented as (N x N) unitary matrices of the Lie group U(iV). 

Likewise, we shall be restricting the configuration space of the plaquette variables (U^ v (x))ijki 
to be compact. In the present case, however, there are four viable definitions of the compactness 
for plaquette variables. The first one, which we call as 'rank-/V 3 class' is defined by 

rank-TV 3 : £ (u^ (*)) . ^ (u^ (*)) ^ = ^irfjfhu 

E(^W)^(^ (x) ) w = /7 5 "'^' 5 "' 



9 




Figure 3: Chan-Paton factor assignment for the cube corresponding to 11 [C„ v x(x)] in Eq. O.41) . 

E(^(x)) yH (E^(x))* /w> = is y/ a*8„,. 

The factor ^ in the right-hand side was introduced by taking an appropriate overall normalization 
of the plaquette variables U fAr (x). All four conditions are necessary, else the discrete rotational 
symmetry of the hypercube lattice would not be warranted. 

For 'rank-Af 2 classes', there are two possible choices. We will label them as N 2 (T) and N 2 (ll) 
classes. They are 

rank-iV 2 (I) : W^vW) (M*)V,, = 8 «%' ( 3 - 8 ) 
k i \ ' ijkl V / i'fkl 

Z{vA*)) m {vAx))] m = W* (3-9) 



'■./ 



ifcvW). J^W)! ( . H) = 8*8//', (3-H) 

TT \ / ijkl \ / i jkl 

Jl K 



and 



rank-^V 2 (II): £fe v (x)) fev(x))* , = 8 /( v8 w 

\ / IJK/ V / 171c! 



(3.12) 
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respectively. 

Finally, there is the 'rank-/V class', defined by 



rank*: E (M*)) (M*))* = 



"jkl 

ifkl 
* 

* 

ijkV 



N8 Jf 



N8 W 



(3.13) 
(3.14) 
(3.15) 
(3.16) 



Notice that each conditions are related by lattice rotational symmetry, based on relation to the 
classical continuum limit. The four classes of compactness conditions are depicted geometrically 
in FigS 




1 



1 



(a) 



(b) 



(c) 



(d) 



Figure 4: Four classes of compactness conditions: (a) rank-N 3 conditions, (b) rank-N 2 (I) condi- 
tions, (c) rank-N 2 (II) conditions, (d) rank-N conditions. Horizontal and vertical directions refer to 
jii and v, respectively. Chan-Paton indices at each overlapping edge are summed over. 



Let us define the configuration space of the plaquette variables obeying these conditions as 
X N 3,X N 2^,X N 2(jiy and X#, respectively. Obviously, 

X^3 C Xjy2(j) and X^^ C X#. 

It thus seems that four inequivalent definition of the theory exists at the quantum level. 

Not all of them seem, however, physically acceptable. A condition we may impose is that the 
plaquette variables with a given compactness condition are always reducible to the link variables 
in a suitably prescribed limit. For example, consider 'dimensional reduction' § of, say, the d-th 
direction with truncating the degrees of freedom of the internal indices as 

{Udn{x))ijki => §ji(U M (x)) ik , ( u /jd(x))ijki => 8,-jt(f/ /i (jc) t ) j h (U dd (x)) ijk i => 8 ik 8j h (3.17) 

'"We will discuss more on lattice dimensional reduction in later subsections. 
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where x — (x 1 ,- • • ,x d ^ 1 ). This is a procedure of the standard dimensional reduction accompanied 
with a truncation of the internal degrees of freedom transforming as adjoint representation of V M (x) . 
Then, {U^(x))ik represents a link variable in Wilson's lattice gauge theory as in Fig. [5]obeying the 
standard gauge transformation rule 

U^x) -)• V{x)U^{x)V{x+jl)\ V(x) e U(JV). 



x+\i 



x+, 



x+d 



~k 



x+\i 
k 



Figure 5: The plaquette variable (Udu{x))ijkl reduces a variable on the link (x,x + p) by truncating 
the degrees of freedom of the indices j, I after the dimensional reduction with respect to the d-th 
direction, 

In the lattice dimensional reduction, the tensor gauge theory action Eq. (l3 . II) is reduced to the 
Wilson's action. We require that the compactness conditions are consistent with Eq. d3.7l) under 
the reduction Eq. (|3.17|) . As demonstrated in Appendix A, it is easy to see that the conditions of 
rank-iV 3 and -Af 2 (II) are ruled out, and that the remaining two, rank-A^ 2 (I) and rank-/V conditions, 
are consistent. 

Henceforth, in this paper, we will consider the two classes for the compactness conditions: 
rank-./V 2 (I) conditions and rank-Af conditions. In the rest of this paper, we consider exclusively 
the rank-A^ 2 (I) class. As will become clearer, this class of theory yields a physically meaningful 
classical continuum limit under the standard dimensional reduction keeping all the internal indices. 
For the rank-Af case, discussions are devoted to Appendix B. It turns out not to lead to a physically 
meaningful continuum classical theory after the dimensional reduction. As such, we conclude 
that consistency with the Yang-Mills theory via the dimensional reduction and the continuum limit 
therein singles out the rank Af 2 (I) class as the only physically meaningful compactness condition. 



12 



The partition function is defined as 

Z=[tDUe- s , ©t/^nn^vW], (3-18) 

J {x} AVV 

where the measure [dU IJV (x)] is induced by the norm: 




and the normalization is taken as / [dU^ (x)] = 1. It is readily seen that the norm is gauge invariant. 
So is the measure defined from it as well. For N = I case, it coincides with the formulation of 
compact U(l) two-form tensor gauge theory on a lattice [|2TT|. 



4 Relation to Dynamical Quantum Yang-Baxter Maps 

In the previous section, we formulated lattice approach of nonabelian tensor gauge theory. The dy- 
namical variables are elementary plaquettes, which carries not only spacetime orientation indices 
but also 'color' indices attached on four links around each elementary plaquette. 

Before proceeding further, in this section, we point out the lattice theory is closely related to 
so-called dynamical quantum Yang-Baxter maps ll28l . To appreciate possible connection most 
transparently, let us first consider ground-states, equivalently, saddle points of the lattice action. 
From Eq. (l3.4l) . we see that the ground-states are characterized by the condition 

(u, v (x)-U vX (x)-U^(x)) = (u^(x + v)-U vk (x + fi)U^(x + %)) , (4.1) 

\ / kqm jur V / kqmjur 

where left- and right-hand sides are abbreviations of 

\ / ijkl \ / Iqun \ / nrmi 

(u Xm (x + v)) (u vX (x + ju)) (U, v (x + i)) (4.2) 

\ / kptq \ / msp j \ / utsr 

respectively. This is a set of cubic algebraic relations among the elementary plaquette variables at 
each lattice sites x and orientations of each cube. We shall make connection with the dynamical 
quantum Yang-Baxter maps progressively, first taking all plaquette variables constant valued and 
then taking account of dependence to the lattice sites and orientations. 

Denote by [x; an elementary cube based at x and oriented along //, v , % directions . Consider 
on the elementary cube a colored open string whose endpoints are fixed at mutually antipodal 
vertices. We can define the holonomy around the elementary cube by lassoing the string around 
the cube once, viz. nonabelian holonomy under parallel transport of the string around a closed 
manifold of S 2 topology. The operation involves six plaquette moves on the cube: three plaquette 
moves plus their conjugate ones. Each set of three moves altogether is a map X ®X ®X — > 
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X ® X ® X , thus carries 6 indices. This then defines the Yang-Baxter map and the curvature is 
measured by a 'deformation' of the Yang-Baxter equation: 

U, v U vX U mX = e ia3H ^ x U vX U, v , 

for each elementary cube [x; p$%] . We immediately notice that the standard quantum Yang-Baxter 
equation corresponds to the equation for flat two-form connection, viz. vanishing parallel transport. 

The way the plaquette variables are ordered is not arbitrary. Rather, it is intrinsically related 
to the way we define the parallel transport in terms of 'lassoing' the string. To see this, start with 
an 'open' string whose two endpoints are held fixed at the two oppositely diagonal vertices. It 
extends over three links, one per each directions of the cube. By a forward 'move', we define 
parallel transport of two adjacent links across the plaquette the two links belong to. We then see 
that succession of the moves ought to be such that the first three moves and the next three moves 
involve are oppositely ordered. In Figj6[ for instance, the string initially along the lower half cuts 
(composed of three links) the cube into two pieces. We first move the string across the lower half 
plaquettes of the cube: after move in the order of \pv], [vA] and \jjk] plaquettes, the string is now 
located at the along the upper cut (composed of three links). 




Figure 6: Quantum Yang-Baxter interpretation of holonomy around the cube. 



The result of these moves is represented by 

We then continue the move through the upper half plaquettes of the cube. In this case, the move 
is again in the order of [vA] and \jJk], but in the transposed orientation in each plaquette. The 
result of these moves is represented by 

zyt r/t 7/f 
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Thus, putting togeter the two set of operations, a complete lassoing of the cube by a string is given 
by 

U, v U vX U^ v ul x U^ (4.3) 

and is precisely the holonomy of the antisymmetric tensor potential. 

So far, we suppressed dependence of the elementary plaquette variables on lattice sites x and 
orientations. Once we reinstate their dependence, the Yang-Baxter maps become enormously com- 
plicated. The dependence on lattice sites x and on orientations of elementary cubes is not arbitrary, 
however. In particular, orientation of elementary plaquettes is closely tied with directional depen- 
dence. Remarkably, we see that these dependences are put together into the so-called dynamical 
quantum Yang-Baxter equations [?] obeying: 

U, v (g-h X )U pX (g)U vX (g-hV) = U vX (g)U^(g + h v )U, v (g). (4.4) 

Here, g is an element of a Lie algebra Q and h is an element of its Cartan subalgebra. Then, making 
the expansion 

U MV = l+hB MV (h) + 0(h 2 ) (4.5) 
for an expansion parameter h, we find that the leading-order contribution is given by 

d^vA + d v£fy + dxB/w ~ [Bfj» , B vX ] - [B yX , B XfJ ] - [B^, B^] = 0, (4.6) 
where the derivative d A , is defined as 




(4.7) 



for the Cartan subalgebra of the Lie algebra Q over C. The left-hand side of Eq. (l4.6l) takes precisely 
the structure we may interpret as the nonabelian tensor field strength. So, we see that the dynamical 
quantum Yang-Baxter equation is the equation for vanishing field strength over the base manifold 
spanned by the Cartan subalgebra of Q . 

5 Parametrization of Plaquette Variables 

To proceed further, we will need to parametrize the plaquette variables (t/ uV (x))^. /t/ that satisfies 
the compactness conditions. In this section, we will construct explicitly a parametrization for the 
plaquette variables defined by the rank-/V 2 (I) compactness conditions Eq. (l3.8H3.1 II) . 
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5.1 Lax pair ansatz 

Among the compactness conditions Eqs. (l3.8l - I3.1 II) . Eqs. (l3.8[ 13.101) imply that the plaquette vari- 
able ([/^(x)) rkl is a \J(N 2 ) matrix with respect to (z, j) and (&,/) pair of indices, referred as row 
and column indices. Likewise, Eqs.( |3.9[ 13.1 II) imply the same but now with respect to (/, i) and 
(j,k) pairs, identified as row and column indices. From these two requirements, we assume that 
(Ufwix)) .j k{ , viewed as (N 2 x N 2 ) matrices, are projected onto a direct product of two matrices 
where the one acts on the indices, i and k, and the other on j and / 0. Therefore, the plaquette 
variables are parametrizable as 

(tf/ivto) := (^vW) ft (^(x)) , (5.1) 

viz. as a direct product of W fW (x) and W fJV (x) belonging to U(/V) x \J(N). 

The two matrices are not independent. The charge-conjugation relation Eq. d3.5l) imposes that 
U(N) ~ U(A^), and that W^(x) = W^xy. Thus, the plaquette variable is parametrizable in terms 
of U(A0 matrices U^\x) = W^(x) as 

(^ W ) ( . w =(^' W ) a (^W).,. (5-2) 

Notice that the diagonal U(l) part of t/y (x) for either /j < v or /j > v are redundant for parametriz- 
ing the plaquette variable. Fortuituously, these extra U(l)'s do not interact with the rest, so they 
will be mod out straightforwardly. The gauge transformation of Eq. (l3.6l) gives rise to gauge trans- 
formation of U$fi (x) as 




(5.3) 



We are thus parametrizing the plaquette variable as a square of J-species of U(AT) link variables, 
Uyfi (x) with n = 1 , ■ ■ • , d, each of which transforms as an adjoint under the gauge transformation 
Vfj(x). Intuitively, the newly introduced link variables in Eq. (l5.2l) are interpretable as the 
variables residing on two 'dual' links of the original plaquette, as depicted in Fig|7l 

Actually, the dual lattice interpretation extends further. In terms of the newly introduced link 
variables Uy(x), the gauge-invariant cube of the plaquette variables U [C^yjx)] m Eq. d3.4b are 
expressible as a triple product of 'dual' plaquette, as depicted in FiglHl and becomes 

u[C^(x)] = <uM[P vX (x)} utofaix)], (5.4) 

II As discussed in Appendix C, there are solutions that keep degrees of freedom of the plaquette variables of O (N 3 ). 
However, these solutions also reduce the U(N) gauge symmetry toU(l) , which is not enough to eliminate all negative 
norm states (ghosts). Therefore, they do not appear to lead to unitary theories. 
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X + V 



I 



X + jl + V 



J 



A 



Figure 7: Parametrization of the plaquette variable ({/^(jc)).^ (flie dashed lines) in terms of 'dual' 

links (the solid lines). The d -dimensional link variables (Uy(x))ik, (Up (x)^)ji can be associated 
with the 'dual' links. 



where the d-species of 'dual' plaquette 11 ^ are 



U^[P yX (x)] = tr f/; ; W^(x+v)^' (x + l)f/r'(x 



(5.5) 



Here, 'tr' refers to the trace operation with respect to (N x N) matrix indices of the dual link 
variables. 

Putting together, the nonabelian tensor gauge theory obeying the compactness conditions Eqs. (l3.8l 
- 13.1 II) is defined by the partition function: 



Ztensor = j[^U] exp(— Stensor) 



where the action is 



(5.6) 



(5.7) 



and the integral measure 

!Dt/ = JJ[U(jc)] where [U(x)] = fl' [dU^ (x)]/Vol (u(l) rf(rf - 1)/2 ) 

/i,V=l 



{-v} 



is given in terms of U (N) -invariant Haar measure [dU^\x)] of the dual link variables. The 
prime (') in the sum in Eq. (|5.9b refers to summing over all nondegenerate plaquette orientations. 

** In the followings, for the Haar measure, we will adopt the normalization convention f[dUy(x)] = 1. 
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A A AAA 

x+v+A, X+M-+V+A 




x x+p. 



Figure 8: The original cube corresponding to 11 [C„ v x(x)] (the dashed line) and the three 'dual' 
plaquettes U [P v x (x)], U ^ [^(x)], U ^ [P^(x)] (the solid lines) for the case /u, v, A, all different. 

Likewise, the prime (') in the product in Eq. (15. 101) refers to integrating over all nondegenrate dual 
link variables. Notice that we eliminated the aforementioned over-counting in U<f (x) for /u < v or 
/li > v by dividing the volume of overall U(l) rf ^ -1 )/ 2 gauge groups. 

5.2 classical continuum limit 

Before proceeding further, we will study classical continuum limit of the action Eq.( |5.9l ) and show 
that the continuum action is manifestly J-dimensional Lorentz invariant. 

At this stage, we will make a slight generalization of the theory we have constructed by in- 
corporating into the theory the variable (i/^uC*) W/ which is associated a collapsed plaquette as in 
Fig.lE 

As will become clear momentarily, adding these variables are imperative in order to obtain a 
Lorentz invariant theory after taking a classical continuum limit. Still, the generalized theory is 
defined by the action Eq. (|5.7| ) and the measure Eq. (|5.8l ), where the sum over the Lorentz indices 
is unconstrained and the volume of overall groups is mod out for d(d + 1) U(l) gauge groups. 
Explicitly, the new action is defined by 

w=-PEEfv^i Re (« M [M*)]o (v) [V*)]« w [^W]-^ 3 ) > (5.9) 
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Figure 9: Degenerate class of plaquettes corresponding to the variable ( U M Jx) 



ijkl 



while the new functional integral measure is defined by 

0U = Yl\U(x)] where [U(x)] = f[ [dU^ (x))/Vol (u(l)' /(rf+1)/2 ) . (5.10) 
W ft v=i 

Notice also that the gauge invariance and the parametrization of plaquette variable into split vari- 
ables are straightforwardly extendible to degenerate variables as well. 

To take continuum limit of the classical action, we expand the dual link variable U^\x) of 
U(AT) gauge group as 



UP (x) = W = 1 + ia 2 A^ (x) - y (a® (x)) + 



(5.11) 



Here, Ay (*) ate (iV x A'') Hermitian matrices, where indices and V transform as Lorentz vector 
indices. Accordingly, the U(A0-invariant plaquette Zl^>[P v x(x)] is expandable as 

2- 



U W (PM) = N + iahrtfMx 



'vA. 



^a 6 tr 



+ 0(a' 



fy\x) = A v Af (x)-A x A v W (x) 



(5.12) 
(5.13) 

where A v denotes the lattice difference operator: A v /(x) = ^(f(x+v) —f(x)). Hence, the classical 
action is expanded as 

W = ^ £ £ [tf 3 tf$> « 2 +^tr (xf + $ {xf + $> (x) 



{x}fi,\,X=l 



2g 7 



^ I N^(x) 2 + 3Nhr(f^(x) 2 ) +0(a 7 ) 



?M^2 



+ 0(a 7 ) 
(5.14) 
(5.15) 



a 2 in the exponent of Eq. (15.1U can be regarded as the area of elementary plaquette, the expansion appears 
natural only for /j ^ v, and not for /j = v for which the corresponding plaquette is degenerate and has a vanishing area. 
However, as we will see, the continuum classical action is Lorentz-invariant only if we take the same parametrization 
Eq. d5.lll l for /j = v as well. 
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in the classical continuum limit: 

a^O and -4 = ^ d = fixed. (5.16) 

g 1 

Here, we decomposed the field strength f$ (x) into diagonal U(l) and traceless parts: 

/v ( ?W + = 0. 

The difference operators in Eq.( |5. 13b should be understood as ordinary derivatives in the contin- 
uum limit. 

For the diagonal U(l) part, we have reproduced the totally antisymmetric 3-form field strength 

flJjW = ~4x -/£ )U(1) ( X ) -f$ m ( x )- Thus > in the iV = 1 case, the result Eq.d5T5l) 

reduces the previously known action of abelian tensor gauge theory. For U(AT), the 3-form field 
strength is constructed as an object carrying six color indices: 

The field-strength is manifestly gauge invariant and has the symmetry properties: 

= = = <5 ' 17) 

meaning that H is antisymmetric under permutations with respect to the three sets of indices 
(fi, i, i'), (v,j,f), Ck,k,kf). Utilizing these properties, the continuum action Eq. (l5. 151) is rewritable 
compactly as: 

S = ^ 1 [d d x £ T %3 (^(x)) 2 , (5.18) 

where 'Tr^' refers to the trace for (N 3 x N 3 ) matrices. It is defined for an generic element fl-ijki'fk 1 
as 

T %3 (^)= £ X ljkijk . (5.19) 

ij,k 

Recall that, in defining the plaquette variables U /M (x), we have included those on degenerate 
plaquette /j = v. Had one considered the lattice theory without them, the naive continuum action 
would be the same as Eq. (l5.15l) except that the summation over /u, v, X is now restricted to the 
cases /j, v, X all different. Then, the first term in Eq. (l5.15l) would still be Lorentz invariant, but the 
second term would not be so because the traceless part f$ (x) is not totally antisymmetric with 
respect to /j, v, X. We have deliberately kept the degenerate plaquette variables so that Lorentz 
invariant continuum action Eq. (15.151) is obtainable. 

The action Eq. (l5.15l) in the classical continuum limit is purely Gaussian. It does not necessarily 
mean that the corresponding quantum theory is free. For instance, as demonstrated in |fT9l[20ll2TI . 
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d = 4 tensor gauge theory of compact U(l) gauge group leads to charge confinement as a result of 
instanton effects. The result is in direct parallel to the situation Il22l of d = 3 compact U(l) gauge 
theory. Such effects are, however, entirely nonperturbative, having directly to do with the topology 
of the configuration space, and the continuum theory defined by the action Eq. (l5.15l) is perturba- 
tively free. Given that the elementary dynamical variables are the link variables (defined on dual 
lattice), it is actually instructive to understand why the continuum theory is noninteracting. Recall 
that the gauge transformation function V^x) in Eq. (l5.3l) is a Lorentz vector and is expandable as 

V,{x) = e iaA ^ = 1 + iaA^x) - l -a 2 (A„ {x)f + ■ ■ ■ , (5.20) 

with A^jc) being (N x N) Hermitian matrix- valued. Thus, in the continuum limit, the gauge 
transformation rule becomes abelian: 

A { ^(x)^A^(x)-d v A^x). 

This conclusion, despite being formulated in terms of the 'dual' link variables, marks the signifi- 
cant departure from the ordinary nonabelian lattice gauge theory. In the continuum theory, we are 
exploring small neighborhood of the identity in the lattice gauge transformation, and information 
of global structure of the group is lost. The global structure can be made visible by introducing a 
cutoff and considering the theory with compact variables, just like the lattice action Eq. (|3.1|) we 
started with. 



6 Dimensional Reduction on the Lattice 

Before proceeding further, we check an important consistency condition. We shall take dimen- 
sional reduction of the J-dimensional lattice tensor gauge theory and show that the theory Eq. (l3.1l) 
is reduced near the continuum limit to the (d — 1) -dimensional lattice gauge theory coupled to an 
adjoint scalar field. 



6.1 dimensional reduction for lattice vector fields 



We first explain what we mean by 'dimensional reduction' in lattice gauge theory. For definiteness, 
we shall consider the J-dimensional Wilson's plaquette action: 

d 



Sym = -P£ f Re (^ 

{x}ftV=l 



fJV\ 



-N 



(6.1) 



fJV\ 



UJx)U v (x+jii)UJx + V^U v (x) 



where the summation of /j, v with the prime (') runs over the region /j^v. The link variables UJx) 
are N x N unitary matrices belonging to the gauge group G, satisfying Eq. d3.7b . It is parametrized 
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in terms of Lie-algebra- valued gauge potential A /J (x) as 

U ll {x)=e iaA ^ x \ (6.2) 

where again a denotes the lattice spacing. 

To consider dimensional reduction from d- to id — 1) -dimensions, we take the d-th lattice 
anisotropic, treat its lattice spacing r much smaller than other lattice spacing a, and take the limit 
r — > 0. Then, J-dimensional hypercubic lattice collapses to id — 1) -dimensional lattice . Denote 
the d-th link variable U d (x) as 

U d (x) = e irA «M, 

while all others UJx) (ji= 1 , ■ • • , d— 1) remain unchanged in Eq. (l6.2l) . Dimensional reduction then 
goes as follows. One picks up the lowest Kaluza-Klein modes (zero modes) with respect to the 
d-th lattice direction, so that x now represents a site of (d — 1) -dimensional lattice (x\, ■ ■ -,Xd-i). 
In the limit r — > 0, d-th links are shrunken to a point, and the variable U d (x) is now associated 
with sites. We take the limit with the combination <&(jt) = ^Ad(x) kept fixed. The scalar field <&(x) 
represents an adjoint scalar field carrying mass-dimensionrn 
The dimensionally reduced lattice action becomes 

5ym = -PI; £' Ro(u[P, v (x)]-n) -2p££Re(«[P^(*)]-tf). 

{x}ftV=l {x}n=l 

After renaming 

U{x) = e ia ^ = U d {x), 

the action in the classical continuum limit reproduces the known result of dimensional reduction: 

Sym = -^U d - l xtv £ \ (i>(*)) 2 + £ \ {D^(x)f 
Sym j L",v=i 4 ^/=i z 

(6.3) 

= a^A v -a v A^ + /[A^,A v ], D M ® = d M <f> + i[A^<t>}. (6.4) 
The limit a — > is taken while holding 1 /g Y M = 2|3a 4 ~( rf ~ 1 ) fixed. 

6.2 dimensional reduction for lattice tensor fields 

By taking the same procedure, we now examine the dimensional reduction of our proposed lattice 
tensor gauge theory. Scaling the d-th dimension differently, we write the elementary plaquette 



^Dividing by a in the definition of <I>(x) is to keep this as a canonical dimension. 
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variables as 

t/uv (x) = 1 + ia B^ (x) H , 

(6.5) 

Udfi(x) = I + iarB dll (x) + ---, 

(6.6) 

Uftd(x) = I + iarB fjd (x) + -~, 

(6.7) 

U dd {x) = I + ir 2 B dd (x) + ---, 

with V = 1 , • • ■ , — 1 and Ijjki = 8^8^ . As r — >• 0, extending the reasoning given above, we see 
that the variables U dfl (x), U^ d (x) become associated to links, while U dd (x) is associated to sites. 
They lead to vector and scalar fields on the (d — 1) -dimensional hyper-cubic lattice, respectively. 
Both fields carry the mass-dimension one, and the limit r — > ought to be taken with newly defined 
fields 

r 2 

^(x) = rB dM (x), Vpix) = rB^ d (x), &(x) = — B dd{x) 

kept finite. 

After the dimensional reduction, the gauge transformation rules become 




where V(x) is a U(/V) matrix-valued field residing at the site x. From the transformation rules, 
we see that (U M (x))ijki is a variable defined on the link (x,x + fi), transforming in adjoint repre- 
sentation ofV/j(x) and carrying additional internal symmetry group labelled by j, I indices. Like- 
wise, we see that (U^x))^ represents a variable defined on the link (x + £i,x), transforming in 
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adjoint representation of VJx) and carrying additional internal symmetry group labelled by i, k 
indices. The variable (U (x))im is defined at the site x, thus represents a scalar field transforming 
as (adjoint) x (adjoint) of the additional internal symmetry group. 

For the theory under consideration, since the plaquette variable is given in split form Eq. (l5.2l) . 
the fields Ud M , U^, Ujj are parametrizable as 



(6.14) 



(6.15) 



u dd{x) ) iju = (^j^^)^^))^^)^. 
and the gauge transformation rules are 

U M (x) -)• V{x)U,{x)V{x + [i)\ 

U^(x) -> V,(x)U^(x)V,(x)\ 
U(x) -»■ V(x)U(x)V(x)\ 



(6.16) 



(6.17) 



In this case, the vector gauge fields with four indices (l/„)ijM, {^^ijkl separate into the purely 
vector field degrees of freedom (A^)^ and the adjoint matter one (AW) ,7. Also, the scalar {Q?)ijki 
transforming as (adjoint) x (adjoint) splits into two cp's. In the classical continuum limit, the gauge 
transformation by VJx) becomes invisible, so the field (A w ) /7 do not transform. Notice that, to 
fix redundant degrees of freedom of the overall U(l)'s in the split form, one may take 

trAW( x ) = trcp(x) = 0. (6.18) 

The dimensionally reduced action is then given by 

S = "Pi I Re[«[C M W]-/V 3 ] 
{x}fi,v,X=l 

(6.19) 

d-\ d-1 

3p£ £ Re [u [C^ d {x)} -N 3 } - 3p£ £ Re [u [C>(x)] -N 3 } , (6.20) 

{x}n,V=l {x}fi=l 

where the term (zi [Cddd{x)\ —N 3 ) is suppressed because it vanishes trivially as a consequence of 
the split form of U dd (x) ■ In the classical continuum limit, the first term becomes Eq. (l5.18l) . but with 
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d replaced by (J— 1). For the second term, 11 [C^^jc)] consists of the three factors U^[P vd (x)], 
U ( v ) [Pduix)] , U ^ [P/jv (x)] . The last factor is nothing but the Wilson's plaquette action giving 



uW[Pp(x)] =N + ia 2 tr[A^A v (x)~A v A^x)} + a 4 tr [—F^ {xf ] + 0(a J 
The field strength is defined as in Eq. (|6.4|) . The first factor leads 



U^[P yd ( x )} = N + a 4 tr 



+ 0(a 5 ), 



where the (a 2 ) -contribution z'a 2 tr(A v A^ (x)) vanishes due to Eq. (l6.18l) . Thus, the second term in 
Eq. (l6.20l) becomes 



3p£ £ Re[u[C, vd (x)}-N 3 } = 3pW£ £ tr 

{x}//,v=l {x}^,v=l 



-f,vW 2 + (a,a( v )w) 



+ O fa- 



Notice that the contribution is of order 0(a 4 ), overwhelming the 0(a 6 ) contribution of the first 
term consisting of the tensor fields alone. The third term is similarly computed to give the kinetic 
term of the adjoint matter cp(x) at (a 4 ). It does not involve nonlinear coupling to the gauge field 
A jU (x) or the scalar field cp(jc), since such couplings are either absent or are higher-orders in the 
continuum limit. Putting them together, we arrive at the classical continuum action 



S = 



3N 2 
Sym 



[d d - l xtr £ {^fcw) 2 +kA (v) (x)) 2 } + "f(z),,(p(x 

H,V=1 ^ J /j=i 



(6.21) 



with g Y M = P^ 5 d an d m e covariant derivative D^cp is defined as in Eq. (l6.4l) . The action Eq. (l6.21l) 
describes the U(A0 gauge theory with adjoint matter, accompanied with (d — 1) copies of free 
decoupled fields. 

The classical continuum action Eq. (l6.21l) is Lorentz invariant. We emphasize that this result 
is far from being obvious. Since the action Eq. (|5.18b emerges from the 0(a 6 ) part of the lattice 
action, we need to keep the expansion Eq. (l5. Ill) up to the (A^) 3 terms. On the other hand, to 
get the action Eq. (|6.21|) of the order 0(a 4 ), we need to expand U M , U up to (A^,) 4 , (A^) 4 , 
(f> , respectively. Because the mass dimension of fields is modified from taking the dimensional 
reduction (which involves taking factors of lattice spacing a to the fields), the latter case needs 
information of one higher-order compared to the former case of Ay . Thus, it is highly nontrivial 
that the continuum limit Eq. (l6.21l) yields Lorentz invariant action Eq. (l5.18l) . 



7 Strong-Coupling Expansion 

An attractive and promising feature of lattice formulation is the feasibility of exploring nonper- 
turbative physics such as dynamical mass generation and confinement etc. One such method is 
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the strong-coupling expansion, which was applied successfully for the Wilson's lattice gauge the- 
ory, and amounts to expansions in powers of the inverse coupling. Strong-coupling expansions 
are defined intrinsically on a lattice and cannot be derived directly for the continuum counterpart. 
As such, one typically supplements the strong-coupling expansions with a suitable methods for 
extrapolating the results to the continuum limit. Nevertheless, even for finite lattice spacing, the 
strong-coupling expansions can lead to new insights by revealing dynamical mechanisms which 
are typical for strongly interacting continuum theories. With such motivation, in this section, we 
develop the method of strong-coupling expansions for the lattice tensor gauge theory. 

Much as in Wilson's lattice gauge theory, the lattice tensor gauge theory admits gauge-invariant, 
nonabelian Wilson surface operators - a direct counterpart of the Wilson loop operators. Correla- 
tion functions involving these Wilson surface operators are the main interest to us. Hence, we shall 
apply the strong-coupling expansion analysis to correlators involving n Wilson surface operators 
for n = 0, 1 , 2, and extract information regarding free energy, internal energy, and surface tension. 

Remarkably, we will find very distinctive behavior of the strong coupling expansion. For or- 
dinary lattice gauge theory, it is well known that the strong-coupling expansion has a finite radius 
of convergence ||29ll . Here, for lattice tensor gauge theory, we will find that the strong-coupling 
expansion is not absolutely convergent but an asymptotic series in the large N limit. As weak- 
coupling perturbation theories give rise to asymptotic series, we conjecture that strong-coupling 
expansion of the lattice tensor gauge theory at large N is dual to an another, weakly coupled lattice 
theory. 



7.1 gauge-fixing 

On a lattice of finite volume and finite spacing, the total number and the domain of integration 
for plaquette variables are finite. Therefore, the functional integral is well defined even without 
gauge fixing. Still, for a given computation of physical quantities, it is often advantageous to fix 
a suitable gauge. We will make use of this gauge-fixing freedom such that some of the plaquette 
variables are set equal to a prescribed value. Recall that we have included in the set of elementary 
plaquette variables those associated with degenerate plaquette, ujf l \x). We found it convenient to 
use the gauge freedom and set them to unity, uffi (x) = In for all x. To show that this procedure is 
always possible, it suffices to find an appropriate gauge transformation function V M (x) such that 

V M (x)U^\x)V M (x + /u) f = I N . 

We find that 



In 



for x^ = 
for Xy > 
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In the gauge choice Ujf\x) = 1, the lattice action Eq. (|5.9l) reduces to 

Stensor = "PL £ ' ^ ( « W «] « M Wl « W ft" «] " ^) , (7.1) 

{*}//, v,A.=l 

where the primed sum over the Lorentz indices runs over /j, v, X all different. 

Notice that the gauge-fixing we have made is partial: the action is still invariant under a class 
of gauge transformations V M (x) that leaves the gauge choice = 1 intact. Such gauge transfor- 
mations are the ones independent of x^ coordinates, since 

V^{x)V ll {x+ l if = 1 viz. V tl (x)=V tl (x+/u). (7.2) 

Notice also that the gauge-fixing does not introduce any nontrivial Jacobian or Faddeev-Popov 
ghosts either. 

7.2 character expansion 

To proceed further for the strong-coupling expansion, we perform the character expansion lf23~Tl . 
The gauge-fixed action Eq. d7.ll) is a triple product of U(N) characters in fundamental representa- 
tion. So, for £/, V, W G U(iV), we expand as 

exp(3P[xn(f/)Xn(V)Xn(W)]+c.c.)= £ C^^x^Wxr^Xr^W), (7.3) 

Ri,R 2 ,Ra 

where %r(U) refers to the character of the irreducible representation R of gauge group U(N). The 
character of the complex conjugate representation R is related to it as % R (£/) = Xr(U) = %r(U^). 
For the trivial representation R = 0, %o(U) = 1, and, for the fundamental representation R = □, 
%□(£/) = till. The sum over Ri,R 2 ,R3 in Eq. fl7.3l) is for all unitary irreducible representations of 
the gauge group \J(N). The expansion Eq. d7.3l) shows that the expansion coefficient Cr 1) r 2i r 3 ((3) is 
totally symmetric under permutations of R 1? R 2 , R3, and that C R1) r 2i r 3 ((3) = Cr^r 2i r 3 (P) . 
Recall that the characters can split or join under the U(N) group integrals as 

J[dU}%R(UAU^B) = 1 Xr (A) X r(5), (7.4) 

/[dt/]% Rl (t/A)x R2 (t/ t 5) = Sr * ^(A*). (7.5) 

Here, cIr = % R (1) denotes the dimension of the representation R. For example, do = 1, da = d n = 
N. The orthogonality relations of characters (the A = B = 1 case in Eq. fl7.5l )), the inverse relation 
of Eq. fl7.3l) is readily obtainable: 

C Ri r 2 r 3 = J [dU] [dV] [dW] Xr, (U) Xr 2 (V) Xr 3 (W) exp (30 \Xn(U)Xn(V)Xn(W)] + c.c.) . (7.6) 
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Factoring out the contribution Cooo(|3), we express Eq. (|7.3| ) in a more convenient form 



exp 3(3 [%n(U)x n (V)% n (W)} + c.c. = C 000 ((3) 



1 + E'Cr 1R2R3 

Rl,R 2 ,R3 



, (7.7) 



where CrjR 2 r 3 ((3) = CrjR 2 r 3 (P)/Cboo(P), and the primed sum runs over all unitary irreducible rep- 
resentations except (Ri,R2,R.3) = (0,0,0). 

We have computed the character expansion coefficient Cooo(P) from Eq. fl7.7b in Appendix D. 
The result is in power-series of [3: 

Qoo(P)=I>!(3|3) 2 "+ t K) 3 n!(3(3) 2 ". (7.8) 

n=0 n=N+l 

It shows that the expansion coefficient grows as n\ up to the 0($ 2N ), and then decreases by the 
factor < v n < 1 . In Appendix D, we computed the first 20 and 50 terms for N = 2 and 3 cases, 
respectively, and concluded from the result that the suppression by (v„) 3 is sufficient to render 
the power-series Eq. (l7.8l) convergent for finite N. It implies that the power-series expansion of 
Cbnn((3) is convergent as well: 

Cb DD (P) := ^37jpy lnC ooo(P) (V.9) 



= (3(3) l+3(3(3) 2 +13(3(3) 4 + 71(3(3) 6 + 461(3p) 8 + 



(7.10) 



For the second equality, we assumed N > 5. 

Interestingly, the large-order behavior of the character expansion coefficients is quite different 
from that of Wilson's lattice gauge theory. For the latter, the character expansion yields 

exp(YXn(t/)+c.c.) = £c R (Y)x R (t/), C B (y) = J [dU]x R (U)ex V (un(U) +c.c). 

R 

The character expansion coefficient Co(y), which is the same as z(j) defined in Eq. (ID.3l) . can be 
computed explicitly. We relegate details to appendix D and quote here the result from Eq. dD.lll) : 



N Jin °° v 

Co(Y)=^(Y)=L^+ I ( 7 -H) 

n=0 ' n=N+l 

Here, the first iV terms coincide with those of , and the power-series converges for any value of 
Y 

Of notable situation is the large-/V limit. For Wilson's lattice gauge theory, as is evident from 
Eq. fl7.1 II) . the first term with N = °° yields a convergent large-order behavior. For the lattice 
tensor gauge theory, however, Eq. (l7.8l) with N = °° is obviously divergent and is not even Borel 
summable. This imparts a significant departure of our lattice tensor gauge theory from Wilson's 
lattice gauge theory. We will dwell on this issue further later in subsection 7.5. 
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7.3 partition function and free energy 



For the gauge ujt = 1, the gauge-fixed partition function is given by 

Z = e -^ 3 d(d-i)(d-2H, z ^ 

Z = [] expfsp^WM^M^aW^Wl + cx.), (7.12) 

J {x}^<v<X 

where refers to the total number of lattice sites, and the measure 2) 'U is the U(/V)-Haar measure 
for the regular plaquette variables (x) (fi^v). Here, we do not consider the volume factor of 
U(l) groups as in Eq. (15.10l) since it merely produces an irrelevant constant factor independent of 
P and N. 

Once the expansion Eq. d7.7l) is made for each cube, Eq. d7.12l) can be written as 

^d{d-l){d-2)9Q 



Cooo(P) 



f^'uu n 



{x} p<v<X 



R] ,R2> R 3 



(7.13) 
(7.14) 



Here, P$(x) denotes the plaquette P v x( x ) formed by the dual link variables carrying superscript 
(ji) , viz. P$ (x) = t/ v W (x) (x+\) C/ V W (x + Xy ufd (x) f . The triple product of the characters in 
Eq. (l7.14l) corresponds to the elementary cube or triple product among its 'dual' plaquettes (depicted 
in Fig. [8]) carrying the representations Ri , R2, R3. The integrals in Eq.( |7. 14b is carried out by making 
use of the integration formulas Eqs. (|7.4[ 17.51) . Nontrivial contributions come from situations that 
elementary cubes, each of which are labelled by the representations Ri,R2,R3 are glued together 
into three-dimensional closed manifolds on the lattice. Consider the simplest one of such cases. 
It is that eight cubes are glued together to form a manifold of S 3 topology. An example of such a 
configuration consists of the following elementary cubes: 



(1): 




(2): 


JCri 


(3): 


%R 2 


(4): 


%R2 


(5): 


XRi 


(6): 


%Rj 


(V): 


%Ri 



to 



P£(x+4) ) %R2 (p[f(x + 4) ) Xr 3 [P£>(x+4) 



,(3), 



^to 



4 2) (*+i: 
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p 

%R3 



(3), 



Xr 3 1^24' W ) 



^ } to 



Pi?(x+t 



,(4), 



Xr 4 [t J i 2 '( x + l . 



4?to 



%R4 



*s?to 



P&\* + ty %r 3 (p[1\x + 2)) %R4 (pit\x + 2; 
/ 5 | 2 ) to)xR 2 (/ 5 1 ( 4 ) to)xR 4 (^i ( 



,(4), 
12 1 



(7.15) 
(7.16) 
(7.17) 
(7.18) 
(7.19) 
(7.20) 
(7.21) 
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(8) : Xr, (P${x+fy % R2 (p|?(x+3)) % R4 (/^(x+3)) , (7.22) 



where each cube is represented in terms of characters of 'dual' plaquettes, and hermitian conjuga- 

1 



tion relations such as P$ (x)^ = (x) are used repeatedly. The group integration then yields 



[dU] (1) x (2) x (3) x (4) x (5) x (6) x (7) x (8) 



(d Rl d R2 d R3 d R4 ) 4 

Intuitively, the result can be understood in terms of the resulting three manifolds. Rewriting con- 
tribution of each representation as d R =d^x d R , the power '2' in the first factor is interpretable 
as the Euler characteristic of § 2 made of six 'dual' plaquettes corresponding to % R , and the power 
'—6' in the second factor is the correct normalization of %r that permits 't Hooft's large-iV power 
counting. Summing over all possible representations of the cubes, the total contribution is given 
by 

\ 4 ^ ^ ^ 2 

J J J J (CriR 2 R3(P) C R2R3R4(P) C RiR4Ri(P) C R4RiR2(P) ) > (7.23) 
a Rl a R2 d Ri d R4 J v J 

where the prime (') of the summation stands for excluding the term (Ri,R,2,R.3,R4) = (0,0,0,0). 
The leading nonzero term comes from the case Ri = R2 = R3 = R4 = □ and its complex conjugate, 
yielding 2(3|3) 8 //V 16 . The next order contribution comes from the case the representations involve 
identity or adjoint. Because Co,o,ad(P), C , a d,ad(P), C ad , a d,ad(P) all start with (3(3) 2 , it gives the 
contribution of o(|3 12 ). Therefore, we find the partition function in strong-coupling expansion as: 

q»(P))**-"^ 

(7.24) 
(7.25) 



1 + I d{ d — \)(d — 2){d — 3)fA6 • ^5 (Cbo (P)) 8 + (P 



12\ 



Here, as explained already, the series-expansion in Cr[R,r 3 (P) represents sum over closed three- 
manifolds. Since terms higher-order in CrjR,r 3 (P) come from contractions among the same vari- 
ables repeated many times, the power series of 3P in CrjR 2 r 3 (P) is interpretable as contributions of 
singular, degenerate closed three-manifolds. 



7.4 Wilson surface observables 
7.4.1 Nonabelian Wilson surfaces 

We begin with a digression regarding nonabelian Wilson surfaces. It is normally considered that, 
for nonabelian 2-form gauge theory, Wilson surfaces are ill-defined. We now show that the no-go 
theorem is evaded for the class of nonabelian tensor gauge theory studied in this paper. 
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Consider the trace of parallel transport around a closed surface E (3E = 0): 

TtU(Z) = HUp. 

Pel. 

In general, the surface E would be self-intersecting. We will call expectation value of such vari- 
ables: 

<W(E) = (Tr£/(E)), (7.26) 

as the Wilson surface observables, and expectation value of their products: 

W (E b - •-,£„) = (TrC/(E 1 )TrC/(E 2 ) ■ • -Trf/(E n )) (7.27) 

as Wilson surface correlators. The Wilson surface observables measure internal energy of the 
system, and the (connected components of) the correlators measure correlation length, mass gap, 
spectrum, etc. 

Wilson surface observables constitute the fundamental basis of the system. Indeed, extending 
the argument of [31], one can assert that every gauge-invariant operator 0, which depends con- 
tinuously on the plaquette variables can be approximated arbitrarily well by a power-series of the 
Wilson surface variables: 

^c(E 1 ,---,E„)Trt/(E 1 )---Trt/(E, ; ). (7.28) 

n>0{£„} 

An important feature of proposed nonabelian Wilson surface observables is that it bypasses 
folklore that such observables are afflicted by ordering ambiguity. For a Wilson surface observable 
of minimal size, viz. the one defined on a cube, we have already shown that there is no ambiguity. 
Could there be any ambiguity when the observables encompass cubes more than one? We will 
now argue that there is no ambiguity by illustrating a few nontrivial cases. The first one involves 
two cubes, as shown in the figure [KB Although there are two possible routes of color move across 
the plaquettes involved, 

We can also illustrate our claim from more sophiscated string move in a given Wilson surface 
observable, as depicted in FigQTJ 

7.4.2 One-point correlator and internal energy 

Having shown that nonabelian Wilson surface operators are well-defined physical observables, 
we now consider the simplest of these operators, taking a rectangular shape: W X (I,J,K) which 
represents a box with the positive orientation, composed by the three edges (x,x + I\), (x,x + 
72), (x,x + K3). It is constructed by tiling faces of the box with the plaquette variables U^. 
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Figure 10: Two alternative routes around two disjoint plaquettes between two string configurations. 




Figure 1 1 : Two alternative routes around an elementary cube between two string configurations. 
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In computation of the vacuum expectation value (W X (I,J,K)) in the strong coupling expansion, 
contributions correspond to various three-dimensional manifolds bounded by the box W X (I,J,K). 
Let us begin with 1, 1)). It is nothing but the internal energy and is computable as 

= , (i -m-2) km) lnz <7 - 29) 

= C nnn (P)+2(J-3)^(c nnn (P)J +0((3 9 ), (7.30) 

where the first term comes from the minimal volume configuration, and the second from elemen- 
tary fluctuations consisting of seven cubes. 

Extending the computation to general Wilson surface (W X (I,J,K)), we find 

UK r ~ , -i 

(7.31) 



(W X (I,J,K)) = (^M) xN i +J+ k 



V ^ 3 J 



i + 2(j-3)/y^ (c ™f )) +o(p 8 ) 



The leading contribution gives the volume-law, indicating that colored strings are confined at 
strong coupling. The hypersurface tension M represents strength of the volume-law and is de- 
fined by the strong-coupling behavior of the one-point correlator for large I,J,K: 



(W x (I,J,K))~exp(-IJKM -(IJ + JK + KI)a-(I + J + K)t 



(7.32) 



with SI and T being some nonuniversal constants. It is evident from the first term that the one-point 
correlator decays with the volume UK times the hypersurface (membrane) tension M . From the 
result Eq. d7.31l) . the hypersurface tension M is extracted as 

* = m | - 3(3(3) 2 - ^(3(3) 4 - { 4 1 + } (3(3) 6 + o (P 8 ) 

(7.33) 

forN>4. 



7.4.3 Two-point correlators and excitation spectrum 



For the connected two-point function (W x (I,J,K)W x+I 2(I,J,K)*) comi , there are two candidates 
giving the leading contribution. The first one is the case that cubes from the action are all used to 
fill the space between W X {I,J,K) and W X+L ^{I,J,K). There are no cubes filling inside each Wilson 
surface. We call it the case (I) (See Fig. [12]). The contribution amounts to 

/~ \ 2(U+JK+K1)L 



(Case (I)) 



N 3 



1 



x 



l + o(p e 



(7.34) 
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x + I 1 + L 4, 



+ L4 




X 



A 



X + I 1 



Figure 12: Schematic view of a configuration of the case (I) for the two-point correlator 
(tU x (I,J,K) c H^ x+L ^(I,J,K)*)c- Elementary cubes from the action are contracted to fill the space 
between { W X {I,J,K) and W X+I 2(I,J,K) indicated by the dashed lines. 

where the 0((3 6 ) contribution is from elementary fluctuations. 

The second possibility is that the elementary cubes first fill up each interior of the two Wilson 
surfaces and leave holes of the minimal size located at y and y + L4, respectively, and then they are 
connected at the holes via a stack of minimal number of elementary cubes (6L cubes). See Fig. [T3l 
The contribution is given by 



draw a similar conclusion. From these considerations, we find that the theory develops a mass gap 
and screening in that the two-point correlator undergoes a 'phase transition' as the separation is 
varied. The situation is rather analogous to the transition taking place for Wilson loop correlators 
in gauge theory, as depicted in Fig {14] 




(7.35) 



where the overall factor UK comes from the sum with respect to the location of y. 



Comparing the power of Cnnn((3) for the two cases, we can see which configuration dominates 
in the strong coupling. For the case / ~ J ~ K, the case (II) dominates if L > ll i.e. the separation 
is larger than the scale of the surfaces. The case (I) dominates if L < ^1. For generic I,J,K, we can 
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Figure 13: Schematic view of a configuration of the case (II) for the two-point correlator 
{W x (I 7 J,K)W x+I ^(IJ : K)*) conn . The two surfaces •W X (I,J,K) and W X+I 2(I,J,K) are connected 
by the shortest tubular stack of elementary cubes (the dashed lines) attaching at the holes y and 
y + Ll 




(a) (b) (c) 

Figure 14: Behavior of the Wilson loop two-point correlator in gauge theory. At small separation 
(a), the first Wilson loop is joined by a minimal surface to the second Wilson loop. At moderate 
separation (b), the minimal surface deflects in the middle. At large separation (c), each Wilson 
loop is filled up by respective minimal surface, joined by a thin tubular column. 
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7.5 Large-A/ Reduction and Asymptotic Behavior 

In lattice gauge theory, it is well known that large N limit exhibit reduction of degrees of free- 
dom, so-called Eguchi-Kawai reduction. In this section, we will find an indication that a similar 
reduction takes place in the large N limit of the lattice tensor gauge theory. 

To observe an indication, consider again the strong coupling expansion of the partition function 
Eq. (l7.25l) . There, all the correction terms in the square bracket are suppressed by some power of 
1 /N as the (Cnnn((3)) 8 term. Under N — > °° limit with P fixed, we can discard the correction terms, 
at least at leading order in l/N expansion. We then see that the partition function Z is reduced to 
the following zero-dimensional, unitary three-matrix model: 

Z (Z MM )^ d( ^ 1)( ^ 2K5 (7.36) 

where 

Zmm = C 000 (P) = /[dt/][d\/][diy] e 3 P( trf/ )( trV/ )( trM/ ) +c - c -. (7.37) 

In this limit, the elementary Wilson surface operator fV x (\, 1,1) alone yields a nonvanishing vac- 
uum expectation value Qxo(P). The connected two-point correlators vanish, and thus the large 
N factorization 

(<0 X {1,1, 1,1)) = (^(1,1,1)) (^(1,1,1)) (7.38) 

holds. 

In the strong coupling expansion, the matrix model partition function Zmm still captures the di- 
vergent behavior of perturbation series. As N — > °°, the n\ growth in the series Cooo(P) continues to 
infinite orders. Therefore, the perturbation series is divergent and not Borel summable. Likewise, 
for the free energy In Zmm? the perturbation series behaves asymptotically as n\ (3|3) 2 ' ! . 

The Borel non-summability of the strong-coupling expansion implies that the entropy in defin- 
ing the lattice partition function of the tensor gauge theory grows much faster than that of the 
Yang-Mills theory. Intuitively, we can understand this by the following geometric considerations. 
At large N limit, the strong coupling expansion of Yang-Mills partition function is interpretable as 
sum over random surfaces. Likewise, at large N limit, we found above that the strong coupling 
expansion of nonabelian tensor partition function is interpretable as sum over random volumes. 
Thus, we would like to compare the entropy of random volumes in comparison with the entropy 
of random surfaces. So, consider the partition function of random surfaces E of area A. For a fixed 
area A, the partition function takes the form 

Z(A,v,N) =£e-^A-*ArX (7.39) 
x 
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for a fixed area A. The sum is over all possible topologies of the random surface E. We are 
interested in the density of states, equivalently, density of states. At fixed N, it is known that the 
partition function scales as ll34ll 



Z-exp(AlogA). (7.40) 

This implies that the entropy S of the random surface for a spherical topology must behave as 

5(§ 2 )~e cA A v (7.41) 

for some numerical factors c,v. that the entropy of random volume of S 3 topology behaves as 

S(E> 3 ) -exp(oAlogA). (7.42) 

The proof goes as follows. Consider the map from M 3 to itself. The map may be represented by a 
vector X = X(^ 1 ,^ 2 ,^ 3 ). The map has folds at critical points where 



det (lf) =a (? ' 43) 

But, the boundaries of the folds are nothing but 2-dimensional random surfaces. For the latter, we 
already know that the partition function behaves as Eq. (17.40b . Hence, we know that the entropy of 
random volume ought to behave the same way, viz. 

Z(§ 3 )~exp(V 3 logV 3 ). (7-44) 

The reduction implies ultra- local nature of the theory, with no degrees of freedom propagating. 
In Wilson's lattice gauge theory case, such a reduction is different from the large N reduction of 
Eguchi-Kawai ll24l . and does not lead an interesting theory. In fact, the free energy simply gives 
lnz(y) = y 2 in the limit N — > °° with y fixed. In the tensor theory case, however some nontriv- 
ial pieces corresponding to singular configurations of three-manifolds seem to remain after the 
reduction. 

Of course, this argument on the large N limit is rather formal. Discarding the correction terms 
needs to be justified by a careful treatment. About this, we will report elsewhere. 



8 Discussions 

We constructed a lattice model for a nonabelian generalization of two-form tensor gauge theory 
starting from an observation of the Wilson's lattice theory. Requirement of consistent dimen- 
sional reduction to lower dimensional Yang-Mills theory uniquely determines the theory, which 
satisfies compactness conditions of rank-/V 2 (I). Also, computation of strong coupling expansion 
was done for the partition function and some Wilson surface observables. There, the character 
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expansion coefficients much more rapidly grow as a power series of the coupling (3 compared to 
the ordinary lattice gauge theory case. It possibly suggests that singular contributions to random 
three-dimensional geometry are much more dominant than those to two-dimensional one. 

There are many interesting and important related to this work. We will mention some directly 
related problems of them. 

• The continuum limit at quantum level needs to be understood better. In the classical contin- 
uum limit, we found that the tensor gauge theory becomes purely Gaussian. However, there 
is a possibility that the lattice theory we formulated may have a nontrivial ultraviolet fixed 
point, where an interacting quantum theory of nonabelian tensor gauge fields can be defined. 
It would be very interesting to explore the possibility via numerical simulation. 

• Universality, which is also related to a suitable continuum limit, need to be understood better. 
Namely, is the continuum theory independent of the latticization? Here, we constructed the 
model on a hypercubic lattice, where the plaquette variables (Uf^(x))ijki encode the square 
lattice structure (and thus the hypercubic one) to the four color indices z, /. Since the 
dynamical variables explicitly depend on information of the latticization, at a glance the uni- 
versality issue seems to be problematic. (In the ordinary lattice gauge theory case, variables 
are assigned on links. Note that each link variable does not explicitly depend on a lattice 
structure — triangular lattice, square lattice, and so on.) By solving the rank N 2 (l) condi- 
tions, however we can rewrite them as 'dual' link variables. Now, the hypercubic structure is 
not explicitly visible any longer in each 'dual' link variable! Thus, the theory with rank iV 2 (I) 
conditions seems hopeful also from the viewpoint of universality. Anyway, check of the uni- 
versality is an important issue related to the Lorentz invariance of the resulting continuum 
theory. 

• The behavior n\ (3(3) 2/? in the character expansion coefficients reminds us of higher order 
behavior of weak coupling perturbation series in quantum field theory, rather than the strong 
coupling. It might suggest an interesting possibility that strong coupling region of the large- 
N theory is dual to some perturbative field theory. 

• We discussed a large N reduction in somewhat speculative way. Related to the duality, it 
would be interesting to investigate weak coupling phase in the unitary three-matrix model 
similar to Ref. 11251 . 

• Supersymmetric extension of the lattice tensor gauge theory is very important, especially, 
in the context of (5 + 1) -dimensional (2,0) theory. In recent years, there has been enormous 
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progress in formulating supersymmetric Yang-Mills theory on a lattice. In the present con- 
text, however, there are further stumbling blocks that need to be overcome. In (2,0) theory, 
degrees of freedom involves not the whole of tensor gauge field but only self-dual part of 
it. Self-dual tensor field in (5+1) dimensions is chiral, described by field equations which is 
first-order in time. Therefore, once put on a lattice naively, the self-dual tensor field would 
faces the problem of species doubling — the tensor field on the lattice would involve not 
only self-dual part but also anti-self-dual part. The situation is exactly the same as chiral 
bosons and chiral fermions in (1+1) dimensions. 
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Appendix 

A Reduction of plaquette variables to link variables 

In this appendix, we show that the nonabelian tensor gauge theory proposed in this paper is re- 
duced consistently to Wilson's lattice gauge theory by truncating internal degrees of freedom, 
which transform as adjoint representation of V^(x), according to the standard dimensional reduc- 
tion Eq. (l3T71) . 

Under the dimensional reduction, the action density of an elementary cube along d-th direction 
is reduced to that of an elementary plaquette: 



«[<W*)] => N 2 u[P^{x)}, 
u[P^{x)} = tr (u M (x)U v (x 
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where we have used the unitarity relation 




Both of 11 [Cfjdd(x)] and 11 [Q^(jc)] reduce to N 3 , but do not contribute to the action. 
Thus, after the lattice dimensional reduction, the action becomes 

d-l d-l 
Educed = -3(3^V 2 £ £' Re(w[^ v (x)]-AM-p£ E Re(«[C M (x)]-iV 3 ). 
W/i,v=l v y {x}fi,x,X=l 

(A.l) 

Here the prime (') refers to sum over /j, v omitting /j = v terms, viz. 11 [P^x]] = N contribution. 
Evidently, the resulting action describes Wilson's lattice gauge theory together with the nonabelian 
lattice tensor gauge theory, both in (d — 1) dimensions. 

As shown in section 3, in the classical continuum limit, the second term in the reduced action 
scales as 0(a 6 ). On the other hand, expanding the link variable U^{x) — e iaA ^ x \ the first term 
produces the Yang-Mills action at 0(a 4 ). Thus, in the classical continuum limit, the first term 
dominates over the second, and the system is reduced to the ordinary U(A0 gauge theory. This 
conclusion is valid for theories defined by both rank-./V 2 (I) and rank-iV compactness conditions. 

B Alternative nonabelian tensor gauge theory 

In this Appendix, we shall consider alternative proposal for the nonabelian tensor gauge theory and 
study its properties. The alternative one is defined in terms of the rank-/V Eqs. (l3 . 1 31 — 13 . 1 61) . We 
shall construct the theory explicitly, work out classical continuum limit and dimensional reduc- 
tion thereof, and demonstrate that this alternative theory does not lead to a physically meaningful 
theory. 

B.l polar parametrization of plaquette variable 

We will first develop an explicit parametrization of the plaquette variables that solves the compact- 
ness conditions Eqs. (l3. 131 - 13. 161) . 

We find it convenient to introduce a set M of (N 2 x A^ 2 ) matrices, and treat the plaquette 
variable (f^/vM)^/ as an element of M where and (kj) are interpreted as column and 
row indices. For general elements A,-^/ and B^i in M, we define a matrix multiplication by 
(AB) ijk i = Zm,n A ijmnB mn ki- Evidently, I ijU = 8 ik 8ji is an identity element of the product. We 
define trace of matrices as TrA = Y^ijMjij- We also introduce a 'twist' matrix Qjki = 6,78^ acting 
on a matrix A^/ as 

(CA)jj k! = A j ikh (AC)jj k i = Aiji k . 
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Using these notations, the condition Eq. (|3.5l) is written as 

CU fP/ (x)*C=U V(1 (x), (B.l) 

where the 'f ' refers to hermitian conjugation for the N 2 x N 2 matrices. 

In the (N 2 x N 2 ) matrix representation, the compactness conditions Eq. d3.13l - I3.16l) are ex- 
pressible as 

EfcWtVW + ). tvI = Nh>, (B.2) 

i \ / IKl K 

k 

jkfx)^) = N8 jf , (B.3) 

ifevW^.vW) , = NS W , (B.4) 

E(^W^W) = N&jf. (B.5) 

Since U lJN {x) is a matrix over C, we parametrize it via (a symmetric version of) the polar decom- 
position: 

U MV (x) = e^^R^x) e^vM 

for every /j, v, where a is a lattice spacing, and H MV (x) is a (N 2 x N 2 ) Hermitian matrix parametriz- 
ing the \J(N ) subset of the configuration space Xjy. R/jv(x) is a positive semi-definite (N 2 x N 2 ) 
Hermitian matrices. The condition Eq. d3.5l) implies that 

CH^ (x) C = -H Vfl (x) , CR/jy (x) C = R Vfl (x) . 

Also, the compactness conditions Eq. (IB.2l — IB.5I) read 

£ ^H^ x ) R ^ e -^( x )\ = £ (e^^R^xfe'^ W ) = A% (B.6) 

In the polar parametrization adopted above, the field R^ (x) is expandable around the identity 

/ as 

R,„ (x) 2 = I- 2a 2 A (UV (x) + a 4 c 2 A /JV (x) 2 + a 6 c 3 A AA , (*) 3 + • • • . (B .7) 

Here, A fN (x) are (N 2 x N 2 ) Hermitian matrices, and C2,C3, ■ • • are real-valued constants. Taking 
the trace of Eq. dB.6l) . we have 

Tr[^ v (x) 2 ] = N 2 . 

At order O (a 2m ) (m: even), it leads 

c m Tv[A tN {x) m ] = 0, 
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and implies that A fJV (x) = when c m ^ 0. Thus, for A pV (x) to be nontrivial, we will need to set 
c m = for m even. 

The coefficients 03,05,- ■■ are not fixable by the conditions Eq. (IB.6l) alone. Here, we will 
consider the simplest case C3 = C5 = • ■ ■ = and discuss the classical continuum limit in the next 
subsection. As will be shown there, the final form of the continuum action does not change even 
when C3,C5, ■ • • are kept nonzero. Thus, the parametrization of U pV (x) is given by 



n v {x) = ei a H ^ x) Jl - 2a 2 A^(x) e$ a H ^ x 



with (N 2 x N 2 ) Hermitian matrices Hfj V (x), A tiV (x) satisfying the constraints 



(B.8) 



(B.9) 



£ VV v ' J lk J k 



^([^(x),-,^!!)^^)]-]) =0 forn>0 

l V v v ' / kikj 



(B.10) 
(B.ll) 



as well as 



CHm{x)C = -H v Jx) 



and 



CAjj V (x)C — A VjU (jc) . 



Due to the constraints, in general, A /JV (x) cannot move independently of H /jV (x) 



(B.12) 



B.2 classical continuum limit 

The plaquette variable U^ v (x) in Eq. (IB.8l) is expanded around the identity / as 



U lN {x)=I +a z 



iH M v(x) -Ap,(x) 



-H^(xf - -A^x) 2 - l - (H /JV (x)A uv (x) +A /jV (x)H^(x)) 



~H^(x) 3 - -A^v(i) 3 - l - (H^A^ixf +A fiV (x) z H fN (x)) 



+1 (HptixfAptix) +A /jV (x)H^(x) 2 + 2H^(x)A^ v (x)H IJV (x)) 

o 



+ 0(a 8 ). 



(B.13) 

(B.14) 

(B.15) 

(B.16) 
(B.17) 
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From the gauge transformation rule Eq. (|3.6l) and the gauge function expanded as in Eq. (15.20|) , we 
obtain infinitesimal gauge transformation rules for H^ v (x) and A fJV (x) 

(h^ (x) J ->• (h^ (x) j + Si*; (a^A v w J - 5 /7 (a v (x) ) a , 

(B.18) 

(a, v w) ;/H (vw) i7H , 

and observe that A fJV (x) is gauge invariant. Again, as for the theory considered in the text, non- 
abelian interactions at finite lattice spacing disappears in the classical continuum limit and the 
gauge transformation rules are reduced to abelian ones. So, in the continuum limit, gauge invariant 
field strength is 

with the same symmetry property of indices as Eq. (l5.17l) . 

Substituting the expansion Eq. (IB.17l) into Eq. (l3.1l) . after some algebra, we arrive at the contin- 
uum action 

Sensor = \^ £ £ { T %3 (x)^ vX (x)] + 3NTv [(d^x))' 

{x}fi,\,X=l 

(B.19) 

+ 16T [A^(x),A v ^(x),A v (x)] } +0(a 7 ). (B.20) 

Here, Tr^' refers to the trace for (N 3 x N 3 ) matrices, which is defined for a generic element 
A-ijktfV as T %3(^) = Lij,k&ijkijk- Also, T [A fN (x),Ayx(x),A^ tl {x)] expresses a trilinear interac- 
tion term, which is defined for (N 2 x N 2 ) matrices X, Y, Z as 

T[X,Y,Z]= ^ XijkiYimjnZnkmi (B.21) 

i,j,k,l,m,n 

with the cyclic property T[X,Y,Z] = T [Y, Z,X] = T [Z,X, Y] . The classical continuum limit is taken 
as in Eq. (l5.16l) . and the result is Lorentz invariant and gauge invariant. 

A remark is in order concerning the remainders in the small lattice spacing expansion. One 
might wonder if higher-order terms may yield nontrivial contributions. Even if we keep the next 
order terms and consider 

a 6 c 3 A MV (x) 3 + a w c 5 A MV (x) 5 + ■■■ (B.22) 

in Eq. (IB.7l) . we arrive at the same continuum action. In this case, in solving the conditions Eq. (IB.6l) . 
Eq. (|B.9l ) remains the same but Eq. dB.l II) is replaced by 

£ ([//,vW,A,vW]) . = £ ([^vW^Wrt = 0, (B.23) 

k \ J ikjk k \ J kikj 
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£c 3 (a3 v (*))^ = -l£ ([^vW, [^W,A,vW]]) i¥ , (B.24) 
£c 3 (a3 v (*))^. = -l£ ([H,vW, [^W,A,vW]]) h/ (B.25) 

In the small lattice spacing expansion, only the first term of Eq. (IB.22l) is relevant, and it leads to 
a 6 Tr(c3A jUV (x) 3 ). Using Eqs. (|B.24[ IB.25I) , however, the contribution vanishes. Therefore, again, 
we arrive at the same action as Eq. (lB.20l) . 

B.3 lattice dimensional reduction 

Following the procedure of section 4.2, we express the dimensionally reduced plaquette variables 
U d/J (x), U^ d {x), U dd (x) in (N 2 x ./V 2 ) matrix notation as 



U M (x) = U d Jx) = e i ^ x \l-2aAJx)e ia 2 H ^ x 



U^x) = U ftd (x) = e i i H ^^/I-2aA fi (x)e i ^ x 

U(x) = U dd {x) = e l ^^I-2ax\{x)e l ^\ 
where the continuum fields of the mass dimension one are defined as 

Hjj (x) = rH dfJ (x) , Afj (x) = rA d/J (x) , 

H M (x) = rH^ d (x) , A M (x) = rA^ d (x) , 

2 2 

3>(x) = —H dd (x) , T\(x) = —A dd (x) , 
a a 

satisfying the conditions corresponding to Eqs. flB.9l - IB.12l) : 

E(^W) fl y t = E(^W) tt y = 0, 

k 



v v ' / ik jk 



(B.26) 



(B.27) 



(B.28) 



(B.29) 



(B.30) 



(B.31) 
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£ ( • • • , [H,{x)Mx)} •••])= for n > 

l Vv v ' J klk J 



(B.32) 

CH lt (x)C=-B fl (x), CA^(x)C = A M (x), (B.33) 



£([*(*),■■■, [*(*),Titt 



II 



ikjk 



£([Tl(x),-,[Ti(x),*(x)]-]) =0 forn>0 

i V v v / / klkj 



(B.34) 



(B.35) 



(B.36) 



C<&{x)C = -*(*) , Ct) (x)C = T| (x) . 
The gauge transformation rules are unchanged from Eqs. (|6.9l — 16. 131) . Introducing the (N x N) 



matrix notation as 



(H,(x)) [jk[ = (h^^+^ix))^ £ff i ? /) (x)=0 ) 



j 

(B.37) 

(A,(x)).. w = (4 il \x)) ik , £A^' ) (x)=trA^(x) = 0, 

it turns out that /z^, transforms as a vector gauge field and Hp ,Aj£ as an adjoint matter field: 

h M (x) -> 7(x)[^(x)-i3 A< ]y(x) t , 

(B.38) 

H C/0 W ^ ^^(x)^, 

(B.39) 

a{/°(x) -» V{x)A { j , \x)V{x)\ (B.40) 

Indices in the superscript (jl) do not transform in the continuum limit. Also, both 4> and r\ trans- 
form as (adjoint) x (adjoint): 

(*«),,„ ^ E, ( v w),( l, «), / (*w) v ^( v ^» t ), t ( l 'W t )„' 

(B.41) 

- i„ ("(iN j? ('»), 7 „( , '(') , ) H N , )„ < b -«) 
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The dimensionally reduced action now reads 



S = "PE E Re[^[C M (x)]-iV 3 ]-3PE E Re[w[C„ vd (x)]-iV 3 ] 
{jc}a/,v,X=1 {x}/j,v=1 



-3p£ £ Re [« [C^(x)] -N 3 ] - p£ Re [tz [C ddd {x)\ -N 3 ] . 



(B.43) 
(B.44) 



The first term gives rise to the same contribution as Eq. (lB.20l) but with the Lorentz indices running 
over 1 , • • • , d — 1 . The second term yields the O (a 4 ) contribution: 



d-\ 



-3p£ £ Re[^[C^W]-^V 3 ] 
{x}ftV=l 



(B.45) 



= pv£ E tr §tfV^to 2 +^E( D M^W 



A'./ 



(B.46) 



-§ E ([^(W J W]rW.C(x)] + [4 w WX B; W]^W 



+2[^(i) I f(i)][f(i)4 ft) (x)] 



+24 £ (a^x))., tr (a^A^x)) } + o(a 5 ), 



(B.47) 

(B.48) 
(B.49) 



i,j,k,l 

where the field strength and the covariant derivatives are defined by 



(B.50) 



D,fly 



V h 



and similarly for D^Ay 7 -*. 



However, contributions from the third and fourth terms start with the o(a 3 ) terms as 



d-l 



-3P££Re[u[C^(*)]-tf 3 ] 



(B.51) 



£ £ 24^a\ mjn (A^ ) ) jl (A [ ; l) ) nm + 1213^(3^^) (Af j ) ;/ (A 
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(B.52) 



+ 0{a 5 ), 



+^{-NTr[(D^) 2 + (D^) 2 } 



-6 



, n /u jntny^n )jl~r\ A fi )nm\^/j ) jl 



(B.53) 



(B.54) 



+6{H < i k) Hf ) +Af ) A^) (4> 2 +r| 2 ^ 



I nqmq 



+3tr (tff > +4 tt W tr (tf <%M +A^>Af 



-3 tr ( ffWtfW) +4») A ^) tr > +A^A^ 



-3tr ( -4*>hM) tr (tf W^*) -a£%M 



(B.55) 



(B.56) 



(B.57) 



(B.58) 



(B.59) 



+ 12 



^mlnjiflu ^)nm~^~^\mlnj{Afj ) nm 



[^VfV[AfU? 9) : 



(B.60) 
(B.61) 



and as 



-p£Re[u[Q^W]-^V 3 ] 

W 



= 8pa 3 £r[r|(*),r|(x),r|(*)] 



(B.62) 



(B.63) 



+ pa 4 £ { - 6 i&rtm&jUk + r\nkminjilk) {<&lqur&wnjq + "(\lqut{(\umjq ~ (<& 2 + if W] 

to 



3 ( < 5 > nA;rai^ > + ^Infon/H y'/ZAr) (^Iqmfomu jq ~\~ ^\lqnu^\mu jq) 



(B.64) 



(B.65) 
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+ 3 (* 2 + 1 fl,n«(* 2 W) 



nqmq 



y&imkffT[lnjm ^\imkn^lnjm) i^jpiq^\kqlp ^\jpiq^kqlp) 



+ 0(a 5 ), 



(B.66) 

(B.67) 
(B.68) 



respectively. Here, the covariant derivative for <I> is defined according to its transformation property 
as (adjoint) x (adjoint) under U(/V): 

( D f&)ijki = ^m + i{h M )^ja + i(h^. f ^ ifkl -i^ im (h^^-i^ki* {h) vr (B.69) 

and similarly for D„T|. In Eqs. (IB.61l - |B.69l) , Einstein convention for Latin indices is assumed for 
simplicity of the notation. 

All the expressions in Eqs. (|B.49[ IB.611 IB.68I) are manifestly gauge invariant except the first 
line in the right-hand- side of Eq. (IB.61l) . where T]-field is acted by the ordinary derivative instead 
of the covariant derivative. To see that it is nevertheless gauge invariant, we will need to consider 
the O (a) -correction contributions to the transformation rules Eq. (lB.40l) . One can solve Eq.( |6.9b 
iteratively with respect to a and get 



h M (x) -»■ V(x)(/v(i)-^)V(i) t +flGi(^(x)) + 



V{x)H!f%)V{x)i + aG 1 {HV%))\ + • • • 



(//' 



(B.70) 



(B.71) 



£ (v,(x)) V(x)4 jl \x)V(xy+aG l (4 jI \x)) 

f,l' 



jf 

2 



*) f ) + 



where the ellipses denote O (a 2 ) -terms and Gi's are 0(a) corrections originating from the V(x) 
transformations: 



Gi(h M (x)) = \\{d,V{x))h,{x)V{x)^ + V{x)h,{x)d,V{x)^ -idJv{x)d,V{x) 



gM j1 \x)) 



Gi(4 J '\x)) 



(d,V(x))H { /\x)V(x)* + V(x)H { /\x)d,V(x) 



{dpV(x))A ( fJ J '\x)V(x) Jf + V(x)4 jl \x)d,V(x 



(;/), 



(B.73) 



(B.74) 



(B.75) 
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For the transformations generated by V^(jc) = e' aA ^ x \ Eq. (|B.72| ) is exact and produces no higher 
order corrections of a. The transformation rules of and T|, Eq. (|B.42| ), are also exact. Under the 
gauge transformation Eq. (IB.72l) with the (a) -corrections retained, one can check that the first 
line in the right-hand side of Eq. (IB.61l ) is gauge invariant (up to irrelevant (a 5 )-terms). 

In the continuum limit, the O (a 3 ) terms in Eq. (IB.61llB.68l) dominate over the Yang-Mills con- 
tribution Eq. (|B.49l ) of O (a 4 ) . The O (a 3 ) terms originate from the trilinear coupling of Eq. (lB.20l) . 
To try to resolve it, if we assumed the mass dimension of r\ being two and rescaled as T] — > ax\ in 
Eqs. (IB.61[ IB.68I) . the contributions would start from o(a 4 ) and i\ijki would become an auxiliary 
field imposing the constraint (Ap )% = 0. In this case, however, the theory is not Lorentz invari- 
ant because of the quartic interactions of Hp that remain in Eq. (IB.61l) . Therefore, the classical 
continuum limit of the dimensionally reduced action does not yield a physically meaningful theory. 

C o (N 3 ) solution for rank-N 2 (I) compactness condition 

Is it possible to find a solution for the rank-/V 2 (I) compactness conditions Eqs. (l3.8H3.1 II) . whose 
degrees of freedom are of order o(N 3 ) at large Af? In this appendix, we demonstrate that such a 
solution can be constructed, and hence demonstrating that Eqs. (13.8H3.1 II) by themselves are not so 
restrictive by themselves. What actually renders the dynamical degrees of freedom reduced further 
to O (N 2 ) is the orientation condition (13.51) . We also demonstrate that, in the naive continuum limit, 
the gauge degrees of freedom is again reduced to (U( 1)) . 

Two of the compactness conditions, Eqs. d3.8l 13.101) . suggests that the plaquete variables are 
parametrizable as (x) = exp[ia 2 T /jV (x)], where T^v(jc) are N 2 x N 2 hermitian matrix- valued ten- 
sor fields. For the moment, for notational simplicity, we shall suppress the indices /j, V and x. Then, 
the hermiticity of the tensor field T is expressed as 



The other two of the compactness conditions, Eqs.( |3.9l 13.1 II) . put further constraints, which we 
shall consider order-by-order in the power-series expansion of the plaquette variables: 



Begin with the condition Eq. (l3.9l) . At the order (a 2 ), it is satisfied by Eq. dC.ll) . The order (a 4 ) 
puts the conditions 



TijM — Tklij- 



(CI) 



1 2 

Uijki = d ik d ji + (ia 2 )T ijk i + - (ia 2 ) (T 2 ) . jU + ■■■. 



j(ik)^j{ki') = q 



(C2) 



k 

Here, each of is a (N x N) matrix- valued tensor field with the notation 
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Then the hermiticity condition, Eq. flC.ll) , is expressed as 



(C.3) 



At next order, 0(a 6 ), no new constraints arise since Eqs. (IC.3[|C.2l) solve automatically up to this 
order. At the order (a 8 ), we get the constraints: 



y / p (im) p [mn) p (nk) p (A:;') (Id') p (im) p (mn) p (ilk) 



2p(mn)p(nk)p(ki')p(im) > ^p (nk) p (kit ') p (im) p (mn) \ q 



(C4) 
(C5) 



We did not find general solutions for Eq. (IC.5l) . Nevertheless, it is easy to see that the following 
two cases solve Eq. (IC.5l) as well as the conditions for the all orders of a 2 : 



( 



(A) : = V 



V 



T N J 



V 



(B): rW=£Vw ( " ) (V- 1 ) 



mki 



(C6) 



(C.7) 



where in the case (A) {t n )* = th , and in the case (B) each of is a (N x N) hermitian matrix. 
The case (A) means that (T^)ji are simultaneously diagonalizable with respect to the indices j, I, 
while the case (B) requires the diagonalizable structure for the indices i, k. 
Conditions that arise at higher order O (a 2K ) take the form 



P\,--,pk-\ 



+ E 



«,m > 1 



n + m = X" 



0. 



(C8) 



n\m\ 



^ j{kqi)j(qiq 2 ) . . . j<(q m ~\i') p(ip\) j(P\P2) . . .j(Pn-\k) 



(C9) 
(CIO) 



It is a straightforward calculation to see that both of the cases (A) and (B) solve Eq. dC.lOl ). Sim- 
ilarly, it is easy to check that the cases (A) and (B) are solutions for Eq. (13 .lib at the same time. 
Note that in the cases (A) and (B), t„ and t\j' have the 0(N 3 ) degrees of freedom respectively. 
(V has the degrees of freedom of sub-leading order o(N 2 ).) 
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The plaquette variables are expressed as 



U ltf = exp[ia L B ltf \, B^ V =Y,B%T A , B%eR, 

A 

where T A are linearly independent basis taking the form Eq. (IC.6l) or Eq. (IC.7l) . The orientation 
condition Eq. (l3.5l) means 

-KM =I<M ■ 

A J A 

(Cll) 

Note that it relates the upper indices k, i to the lower indices /, j for (t^ a J ^ . For the case (A), it 
reads 

— ^Bpv^Vlm(V l )mjtm l A = ^B A fl ^V] m (V X ) m itm^ A 
Am Am 

(C.12) 

which means the diagonal structure also for the indices k, i: 

m 

(C.13) 

with t A m = t A n . The degrees of freedom reduce to 0(N 2 ). Similarly, for the case (B), it requires 
the diagonal structure for the lower indices /, j reducing the degrees of freedom to (N 2 ). 

Finally, because the gauge transformation rule Eq. d3.6l) ought to keep the diagonal structure in 
the cases (A) and (B), we should consider the gauge group (U(1)) JV , a subgroup of U(N). Thus, 
the solutions do not lead to an interesting nonabelian tensor theory. Also, since the gauge degrees 
of freedom do not have enough degrees of freedom to kill the modes of wrong-sign kinetic term 
(ghosts), the solutions do not seem to lead physically meaningful continuum theory at least in the 
classical level. 



D Character expansion coefficient Qoo ( P ) 
D.l evaluation of generating function 

Here, we compute Cooo(P) we needed in the text for obtaining results for the partition function and 
the correlators among Wilson surface operator. Using Eq. (|7.6l ), we can express Cqoo(P) as 



Cooo(P) - £ 



(3(3) 



In 



to ("0 



n2 



'fm( 



tiUtiU 



(D.l) 
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To compute the U(N) group integral in the right-hand- side, consider the following generating func- 
tion: 



z(y) = f[dU]e^ u+u ^ 

y2n 



= I 



to ("0 



\\2 



[dU] trC/trt/' 



(D.2) 
(D.3) 



The large-Af behavior of this integral, with the 't Hooft coupling X = N/y fixed, was investigated 
in [|25l . Here, we will extend it to finite N. Making use of the character expansion: 

e T*V = £ Cr(y) Xr{U) and e y*ut = £ ^(y) %R (t/t) (D . 4 ) 

R R 

and the character orthogonality relations, we find that z(j) — Er ^r(y) 2 - We shall first evaluate 
Cr(y) by carrying out the integral 



(D.5) 



Any unitary irreducible representation is labelled by the ordered integers mi > mi > • • • > m^, 
running over all integer values while keeping the order. Note, however, that the representations 
that appear in Eq. (ID.4l) with nonzero £r(y) are restricted to nonnegative m,'s, i.e. m\>mi>---> 
mpj > 0. When U is diagonalized in the form 

\ 

a f , a g su(^v) 



u = n 



\ 



7 



(D.6) 



in terms of the shifted weights hi=N — i + m,-, the character %r(U) is expressed as 

det J)k (e ih J^) 



Xr(U) 



where 



A(4>) = det M (e 



(D.7) 



Also, the measure becomes 



[dU] 



m 



(D.8) 
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with the normalization for the angular integral as f d£l = 1. Substituting these into Eq.( |D.5l ), after 
some algebra, we arrive at 

n ^ 

j<i 



«v) - nr<*« IJjr -^, - Ilr „, .; ,„ : II (i 



mj + l-j 



(D.9) 



and hence 

z w= n t 2 "" ; n2 - ( ^ nf 1 — 5— V- (d.io) 

m!>m 2 > >m w >0 ( m l ! ) 2 ' • • KM) 2 y V mj + l-jj 

D.2 asymptotic behavior 

By computing the generating function z(y) explicitly up to 0(y 10 ), we were convinced that it is 
expressible as 

N Jin °° 

^) = Itr + I ^ (D11) 

where < v n < 1. In other words, the first Af terms coincide with those of , but remaining 
higher terms are suppressed compared to the corresponding terms of et '. At o(j 2n ), each term 
contributing to Eq. flD.10l) corresponds to a way of partition of the number n into N. For n <N, the 
number of such terms is just the number of partition of n, but as n goes over N it becomes smaller 
than the partition number. This is a reason why the suppression factor arises. In fact, the first few 
terms of v n takes the form: 



1 



ViV+l 



VN+2 



VN+3 



(N+l)\ 

N+l 



(D.12) 



(N + 2)N\ (AH- 2)!' 

(D.13) 

\(N+l)(N + 2) 1 N(N + 3) N + 2 

4~{N + 3)Nl 4(N+l)(N + 2)(N-l)l ~ (AH-3)(AH- 1) ! 

(D.14) 

1 



(AH-3)! 

(D.15) 

As increasing the order, v n gets smaller. In fact, Eq. dD.l II) is consistent with the large-Af result for 
the strong coupling phase studied in E51 



(D.16) 



lim hxz(N/X) = -r-~ for X > 2. 
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From the power-series Eq. flD.l II ), we find that 




for n < N 



for n>N+\, 



(D.17) 



and obtain the result 



N 



Cooo(P) = £n!(3(3) 2 "+ £ 



3 n!(3p) 



(D.18) 



n=0 n=N+l 



Due to the cubic power of the U -integral term in Eq. dD.ll) . the expansion coefficient grows as n! up 
to O (|3 2 ' ! ) . It is quite different from the ordinary lattice gauge theory case, where the corresponding 
quantity in z(y) has behavior of \/n\. 

D.3 Convergence of Cooo (P) 

To understand the large-order behavior of the second sum in Eq. flD.l 81) for finite N, we examined 
the first 20 and 50 terms in the cases N = 2 and 3, respectively. The result is tabulated in Table CD 
For N = 2, the value saturates at n = 5 and then decreases. Suppression rate becomes stronger 
with increasing order. For N = 3, the behavior is similar, and saturates at n = 20. Passing the 
saturation point, the value decrease rapidly. For both cases, the asymptotic behavior shows clearly 
that that series converges. However, the convergence should be considered typical only for finite 
N. As N gets larger, the saturation point shifts to a larger n quickly. For infinite N, the series is 
dominated by the first sum, which is neither convergent nor Borel summable. 
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